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ABSTRACT 


A key element of the U.S. Navy’s transition to an eleetric naval force is an 
Integrated Power System (IPS) that provides continuity of service to vital systems despite 
combat damage. In order to meet subsequent survivability standards under a reduced 
manning constraint, the IPS system must include a fault tolerant control scheme, capable 
of achieving automated graceful degradation despite major disruptions involving 
cascading failures. Toward this objective, online model-based residual generation 
techniques are proposed, which identify explicitly defined faults within a stochastic DC 
Zonal Electrical Distribution System (DC ZEDS). Two novel polynomial approaches to 
the design of unknown input observers (UIO) are developed to estimate the partial state 
and, under certain conditions, the unknown input. These methods are shown to apply to a 
larger class of systems compared to standard projection based approaches where the UIO 
rank condition is not satisfied. It is shown that the partial-state estimate is sufficient to 
the computation of residuals for fault diagnosis, even in such cases where full-state 
estimation is not possible. In order to reduce the complexity of the system, a modular 
approach to Eault Detection and Isolation (EDI) is presented. Here, the innovations 
generated from a bank of Kalman filters (some of them UIOs) act as a structured residual 
set for the stochastic DC ZEDS subsystem modules and are shown to detect and isolate 
various classes of faults. Certain mathematical models are also shown to effectively 
identify input/output consistency of systems in explicitly defined fault conditions. 
Numerical simulation results are based on the well-documented Office of Naval Research 
Control Challenge benchmark system, which represents a prototypical U.S. Navy 
shipboard IPS power distribution system. 
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EXECUTIVE SUMMARY 


The incorporation of an Integrated Power System (IPS) in the DDG 1000 
Destroyer represents a first step in a transformation to an Electric Naval Force [72]. 
Moreover, the U.S. Navy is presently assessing alternative power and control system 
architectures for the Next Generation Integrated Power System (NGIPS) [395]. These 
systems must be survivable and increasingly capable of brokering propulsion and pulsed 
power demands while providing ships service power continuity to vital mission systems 
[91]. In an era of asymmetric warfare and reduced manning constraints, the NGIPS must 
provide automated continuity of service control to vital systems despite major disruptions 
involving cascading failures. Toward this objective, online physics-based fault detection 
and isolation (EDI) techniques are proposed and are expected to operate within a fault 
tolerant control (FTC) framework. The methods seek to identify explicitly defined faults 
within a DC Zonal Electrical Distribution System (DC ZEDS) portion of an IPS. Here, 
mathematical models, together with sensors providing real-time data, are used to validate 
the condition of the system. The results of this research demonstrate various methods of 
generating fault-indicating signatures for classes of sensor, actuator, and plant faults. In 
the case of system actuator EDI, several estimation techniques based on a polynomial 
realization are developed. Numerical simulation results are based on the well- 
documented Office of Naval Research Control Challenge benchmark system, which 
represents a prototypical U.S. Navy shipboard IPS power distribution system. 
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I. 


INTRODUCTION 


In keeping with the vision set forth by the Chief of Naval Operations, Naval 
warships in the Twenty-First Century must have the ability to “ehange, adapt and 
transform” to meet new threats. Indeed, Sea Power 21 defines a Navy capable of 
constantly improving its ability to “project offensive power, defensive assurance, and 
operational independence around the globe” [111]. To achieve the vision, the CNO 
outlines three fundamental pillars; among them is a continuing transformation of 
improved technology. 

Achieving these goals require warships with a higher degree of survivability to 
enable decisiveness, sustainability, responsiveness and agility, which are the four 
fundamental qualities in Sea Power 21. From an operational point of view, survivability 
is the ability to avoid and withstand a hostile environment. From a system-level point of 
view, it can be thought of as a “fight hurt” or “fight-through” capability - that is, the 
warship’s ability to maintain continuity of service despite major disruptions and 
cascading failures in the presence of catastrophic events. 

Within this framework, an automated fault-tolerant control scheme is essential to 
a warship’s fight-through capability. Of particular importance in achieving a graceful 
degradation of warship systems under these conditions is a design strategy that integrates 
combat systems with advanced hull, mechanical and electrical capabilities. According to 
the U.S. Navy Survivability Design Handbook for Surface Ships, developing affordable 
survivability measures under a distributed, dynamically interdependent system of systems 
framework will “constitute the hallmark of the Navy's new approach to reliably automate 
damage control functions while operating with fewer crewmembers” [343]. 

Take for example an Integrated Fight Through Power (IFTP) system, which seeks 
to maintain power continuity resulting from a calamitous exogenous event. According to 
numerous Naval Research Advisory Committee (NRAC) reports and naval power 
systems and automation experts, there is a growing need to develop cost-effective ways 
to solve power automation issues by taking into account the following aspects: 
[390] [391] [91] [264] [72] [76] 
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1. Automation will be essential to ensure survivable and effective surface 
combatants in the warfare environment [271]. 

2. Personnel costs comprise over 50% of operating and supporting costs, 
therefore it is necessary to reduce the number of shipboard war fighters. 
Smart Ship demonstrated that technology and procedural changes can 
improve ship operations [272]. 

3. Well-trained crew members are decision-makers that require real-time data to 
make informed decisions. No longer able to rely on a well-staffed crew for 
data and unable to react to rapid and complex nonlinear system 
interdependencies, an automation strategy must instead respond to cascading 
system faults and failures [392]. 

4. Life Cycle Costs (LCC), which include research and development, can be 
reduced by implementing simulation-based design environments - especially 
in cases where revolutionary technologies are being developed - where no 
historical cost data exists [273]. 

5. Advanced distributed simulations and simulation-based 
design/manufacturing to assess new technologies and ship design options is 
encouraged [273] [392]. 

Given this trend, it is most desirable and advantageous to employ an on-line 
analytically redundant, robust Fault Tolerant Control (FTC) scheme designed to detect 
unpermitted deviations from a system’s expected, acceptable or standard condition and 
perform actions as necessary to maintain continuity of service - in this case, power 
distribution. Herein, the first of a two-step FTC process is presented using robust model- 
based Fault Detection and Isolation (FDI) methods. Characteristics of faults such as 
time, location, size, and relative influence are determined through residual generation and 
analysis methods developed herein. The methods are applied to a DC Zonal Electrical 
Distribution System (DC ZEDS) testbed designed to resemble dc power distribution 
systems currently being researched for installation on future all-electric warships for 
reasons to be discussed. 
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A. ROADMAP TO AN ELECTRIC NAVAL FORCE 

The NRAC report Roadmap to an Electric Naval Force identified several major 
motivations for the navy to shift away from a separated power and propulsion platform 
and towards an all-electric warship [271] [390] [391] [392]: 


• Electric weapons and advanced sensors for superior firepower range and 
resolution require far greater pulse power than currently available. 

• Electric propulsion and auxiliaries possess superior mobility, stealth and 
endurance. 

• Integrated Power System (IPS) allows provision for common electric 
power system for real-time power allocation, reconfigurability, and 
superior survivability. 

• Electric ship platform provides support for off board weapons and sensors 
for superior reach and warfighter sustainment. 


These reasons among many others make transitioning to an all-electric drive 
warship as significant a change as shifting from steam to gas turbine and possibly even 
from sail to steam. With this change comes the challenge of managing far greater power 
requirements - on the order of tens of megawatts, which is only the beginning 
[11][72][90]. Already committed to fielding an electric warship, the U.S. Navy is 
currently building DDG 1000 as the baseline for future DC zonal electric power systems 
[271]. Again, survivability requirements and current trends in power electronics 
technology obviate the need for a power automation scheme capable of providing 
continuity of service under temporal and spatial bursts of natural or hostile disruption 
[390]. A Fault Tolerant Control scheme, then, provides a framework upon which to build 
a desirable, realizable, and tractable solution. 
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B. FAULT TOLERANT CONTROL 


A fault can be described as any deviation from nominal system behavior that can 
negatively impact overall desired system performance. A mathematical model-based 
approach to the detection of a fault, therefore, is sought. From a system automation point 
of view, a representative Fault Tolerant Control architecture can be depicted as in Figure 
I.l [44]. 


/V 


/ 



Figure I.l: Fault Management Architecture 


Notice the fault management architecture contains an execution level, which 
resembles common model-based state feedback control system architecture and a 
supervisory level, which depicts two major processes - diagnosis and controller re¬ 
design. The plant includes modeling uncertainties defined as a vector of unknown inputs 
d . A vector of faults / , cause the plant to deviate from a nominal condition. 
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Assuming linearity, a plant model (depicted in Figure LI) can be expressed in 
continuous time as 

X = Ax + Bu + Ed + Gf 

( 0 . 1 ) 

y = Cx + Fd + Hf 

where state vector xeM", input mg M™, and output ye M'" and matrices A, B, C, E, 
E , G and H are of appropriate dimensions. The terms Ed and Gf model uncertain 
(unknown) inputs and faults, respectively, in the actuator(s) and plant model. In the 
output equation, Ed and Elf represent uncertain (unknown) inputs and faults in the 
sensors. Note that matrices A, B , and C represent the nominal system where all faults 
and uncertainties are assumed to be contained within / and d . It is further assumed that 
matrices E, E, G and H are constant and known, while / and d are typically 
unknown. 

Looking at the reconfiguration level in Figure 1.1, the main objective of the on¬ 
line Fault Diagnosis process block is to determine if the system is subjected to a fault and 
its time of occurrence. Then, all subsequent actions should focus on characterizing the 
fault(s) to the greatest extent possible. These actions should include estimating the fault’s 
location, identification, and severity. If a history of fault activity is needed for 
conditioned-based maintenance or scenario reconstruction, then developing fault statistics 
would be advantageous. This step is represented as the output to the Fault Diagnosis 
process block and is annotated by / . 

The Controller/Plant adjustment process block then receives the fault information 
and performs appropriate controlling action as depicted in the figure. This action might 
be to adjust the parameters of the controller or reconfigure the controller altogether. If a 
fault becomes a failure in the nominal plant, then it might be necessary to operate in a 
reduced capability, damage control reconfiguration mode in order to maintain continuity 
of service. If this is the case, then it might be advantageous to have the Controller/Plant 
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adjustment process block initiate a reconfiguring of the nominal plant and associated 
mathematical models in the Fault Diagnosis block. 

Some work has been done at the execution level within a Fault Tolerant Control 
framework. To a limited degree, fault tolerance can be achieved under a robust or 
adaptive control design. In Robust Control, passive fault tolerance is achieved by a fixed 
controller that is designed to satisfy pre-defined performance specifications. Adaptive 
control methods, on the other hand, are considered active fault tolerance measures since 
controller parameters are adjusted to changing plant parameters (possibly caused by 
faults) as a function of time. In both cases, however, the classes of faults are restricted in 
the sense that they must be “well-behaved” according to their fault definitions and 
mathematical formulations. Adaptive control, for example, performs well only for linear 
systems with slowly changing parameters [44]. Neither method attempts to diagnose a 
broader range of faults nor do they attempt to change the controller or plant structure to 
mask, compensate, correct, re-configure, re-design or otherwise “tolerate” faults. This 
work looks beyond the restricted class of faults to generate residuals at the supervisory 
level. 


C. OVERVIEW AND CONTRIBUTIONS 

As described in Figure I.l, the flow in the reconfiguration level of a representative 
FTC architecture begins with Fault Diagnosis. It is the intent of this work to provide the 
output of the Fault Diagnosis process block - that is, the necessary fault information for 
further processing to obtain the Navy’s goal: achieve automated fight-through capability 
despite major disruptions and cascading failures. Common error masking and error 
recovery methods are poorly suited for a system description under this framework 
[390] [392]. 

The goal of this work is to develop Fault Diagnosis strategies applied to the DC 
Zonal Electrical Distribution System (DC ZEDS) by generating model-based residuals 
from estimates of the partial state. The following chapters will establish a methodology 
for establishing the DC ZEDS models, estimation methods, and Eault Detection and 
Isolation (EDI) methods. 
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Chapter II details the modules of the DC ZEDS system to include how some of 
the nonlinear dynamics are averaged. Since a deductive approach will be taken to 
identify the effects of a system fault, a close observation of fault effects will be 
monitored. It will be shown that output sensor measurements are selected to guarantee 
observability and are consistent with algebraic relationships and conservation of energy 
laws. Chapter III then explains the methodology for choosing estimation methods, 
including several novel polynomial unknown input observer methods. An explanation 
for using linear methods on nonlinear models will also be addressed. Chapter IV 
identifies common mode failures and identifies classes of faults through simulation, then 
provides a methodology for choosing the Fault Diagnosis schemes to developing fault 
indications in the residuals. Finally, Chapter V provides a summary of results. 

Overall, the contributions of this work can be categorized into 2 main areas. The 
first is in a novel polynomial approach to the design of unknown input observers (UIO) 
with particular attention paid to partial-state and input estimation. The second area is in 
the FDI techniques applied to the DC ZEDS application. 

A structural residual set - predicated upon physics-based mathematical model 
estimation - is developed to decouple the effects from a defined set of faults in order to 
deter mi ne fault indications. The unknown input observer (UIO) developed in this 
research will be part of this task. 

Even though UIOs are used for Fault Detection and Isolation (FDI), little 
application research has been done to take advantage of input estimation, which is shown 
to enhance analytical redundancy. In addition, it will be shown that the proposed 
polynomial based approach can be applied to a more general class of systems, compared 
to the standard projection based approach of [175]. 

A polynomial approach to a general UIO for state and input estimation for both 
discrete and continuous time systems is presented for the purpose of fault detection and 
isolation. Unlike UIO methods that are based on a subspace projection of the whole state 
onto a subspace of known state components (such as [73], [80] and [175]), this robust 
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estimation method is based on an estimate of the partial-state representation of the 
system. This is sufficient to the computation of a residual, even in cases when full-state 
estimation is not possible. 
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II. DC ZONAL ELECTRICAL DISTRIBUTION SYSTEM 


Among all shipboard distributed systems, it can be argued that the power 
distribution system is the most vital. If not for its tightly coupled interdependence with 
every shipboard distributed system (including propulsion), it would be in terms of 
potential vulnerability. Indeed, an all-electric drive warship depends absolutely upon its 
power system to provide quality and continuity of service while operating in a hostile 
environment. And with future warships potentially requiring eight to ten times the power 
generated on today’s warships, it will only increase its stock. 

The prototypical power electronics based dc zonal electrical distribution system is 
uniquely designed as an efficient high power density, fast transient response, tightly 
coupled finite inertia system. The sophisticated and robust power converter control 
elements provide nearly-ideal load regulation and transient performance (tens of 
milliseconds) with a high degree of fault tolerance. The system is also characterized as 
having a high degree of automation and nearly instantaneous reconfiguration capabilities. 
For example, a dc/dc buck converter - which is the Ship’s Service Converter Module 
(SSCM) in the representative DC ZEDS system herein - is capable of maintaining a 
constant output voltage regardless of input disturbances. From the output side, this 
attribute is highly desirable for power quality-sensitive loads. 

A. ONR INTEGRATED POWER SYSTEM TESTBED 

The nonproprietary, well-characterized, reduced-scale ONR IPS Testbed DC 
ZEDS model serves as an ideal resource for initial validation of electric warship power 
distribution analysis methods and control schemes. Through sponsorships from the 
National Science Foundation (NSF) and the Office of Naval Research (ONR), the 
MATLAB Simulink based model was developed by a collaboration of academic 
communities from Purdue University and the U.S. Naval Academy. The distribution 
system topology is also representative of the Naval Combat Survivability Testbed, a 
reduced-scale hardware testbed designed to achieve the same purpose. 
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For the purpose of this research, the DC ZEDS model - in both average value and 
detailed model forms - is used to apply novel estimation and residual generation 
techniques for robust model-based fault diagnosis. When considering the topology, the 
DC ZEDS system is spatially designed to improve survivability and redundancy by 
splitting the DC power distribution into a two-bus system - one on the port side, the other 
on the starboard - each being fed by a dedicated power supply. The power distribution 
system is subdivided into 3 zones (numbered forward to aft) that are separated by water¬ 
tight compartments as shown in Eigure II. 1. 



Figure II.l: DC Zonal Electrical Distribution System Topology 


Within each zone, two Ship’s Service Converter Modules (SSCM) are fed from 
the port bus and starboard bus, respectively. Both of their regulated voltage outputs are 
compared and all the zonal load current demand is supplied from the bus with the higher 
voltage. The converter modules have the capability to share power requirements in droop 
mode, but is not considered in this study. Opposing diodes on the lines leading to the 
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load prevent one bus from providing power to the other. The load in each zone is 
identical, consisting of a Ship’s Service Inverter Module (SSIM) feeding an ac load bank. 

If either power supply fails, the remaining bus will supply the entire load without 
interruption. Current limits on the SSCMs alleviate faults at the oring (diode) function 
while some faults within the modules are mitigated through robust controls. It should be 
noted here that the control scheme improves the robustness of the system, but this can 
have a masking effect on detecting and isolating faults. Further discussion on this will be 
addressed in a later chapter. 

As stated before, there are two time-domain simulation models considered in this 
work: the Average Value Model (AVM) and the more accurate Detailed Model. Power 
electronics modeling experts have noted that a so-called ‘detailed’ model is unfortunate 
because the term is rather arbitrary, but have also documented its sufficient accuracy as a 
truth model for initial system analysis and control architecture simulations. For the 
purpose of this dissertation, which is consistent with the literature, ‘detailed’ refers to a 
simulation in which semiconductor switching action is included, even if only on an on or 
off basis. By comparison, the nonlinear average value model refers to simulations where 
the switching is represented on an average value basis resulting in state variables constant 
in the steady state. Although the average value model removes the switching dynamics 
present in the detailed model, both models are comparable in every other respect, which 
includes parameter values, control architecture, dc steady state, and transient response. 

By comparison, the detailed model is the more representative simulation of the 
DC ZEDS system, but the average value model provides the basis for the state space 
models used for Fault Diagnosis. The remaining part of the chapter describes these 
modules in more detail starting with the detailed model, followed by the nonlinear state 
space models of the average value model. It should be noted here that the red markings 
on each circuit diagram represent sensor measurements. By convention, a red sensor 
marking with a negative slope indicates a current sensor, while a positive slope indicates 
a voltage sensor (always measured with respect to ground). Finally, a frequency analysis 
of the two models is presented and a method is proposed to generate residuals based on 
the AVM state space models and measurement data from the Detailed Model. 
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B. MODULES 


1. Power Supply 

a. Detailed Model 

Each 15 kW power supply is fed from a dedicated 3-phase ac source that 
may vary from 480 - 560 V line-to-line rms at 60 Hz and provides 500 V dc output to the 
bus. The power supply consists of an isolating transformer (modeled as part of the 
leakage inductance), a controlled 3-phase bridge rectifier, and dc link capacitance, 
inductance and resistance as shown in Figure 11.2. 




For the simulations herein, the uncontrolled rectifier mode is considered 
with a zero firing angle (a = 0). The output voltage is regulated by a buck 

converter/controller scheme connected to the output of the rectifier. Parameters for the 
system are shown in Table 11.1. 
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Parameter 

Description 

Value 

Units 

h 

Source Commutating Inductance 

1.0 

mH 

^dc 

Rectifier Line resistance (dc link resistance) 

0.03 

Q 

^dc 

Rectifier line inductance (dc link inductance) 

0 

H 

Cr 

Rectifier output (dc capacitor) 

500.0 

//F 


Input ac voltage frequency 

377 

rad / sec 


Buck controller line resistance 

0.03 

Q 

^out 

Buck controller line inductance 

3.0 

mH 

^out 

Buck controller output capacitance 

500.0 

//F 


Table II.l: Power Supply Parameters 


b. Average Value Model 

The average value model for the power supply is developed in [220] and 
makes the assumptions that the rms amplitude of the ac source voltages E is constant, 
only one commutation occurs at a time, and the dc load current is constant (or at least 
near-constant). The average value model makes several reference frame transformations 
of the input voltage to obtain the time-dependent peak voltage of the input, A 

more detailed derivation of the transformation is shown in Appendix B. Since the output 
waveform of the bridge rectifier circuit is periodic over a ;r/3 interval of the angular 
position of the source voltage 6 ^, the average dc voltage can be expressed as 


V 


d 



( 0 . 2 ) 
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With further substitution and simplification, the average dc voltage from 
the rectifier can be expressed as 


V 


d 


3V3 


V2£ 


cos a 


(0.3) 


and the average dc current is a state in the average value model as shown in Figure 
11.3. 




7T 


out 



Figure II.3: Power Supply Average Value Model 

From the model in Figure II.3, the power supply rectifier and buck 
converter are dynamically coupled through a current and voltage source dependent on the 
duty cycle d , which is the control input used to regulate the output voltage . 
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c. Power Supply Control 

The power supply controller design is shown in Figure 11.4 with its 
associated parameters listed in Table II.2. 
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NSC 


{ 4 , 


tf. 


^NSC^ ^ 


(^ Y' 






out 


out ■ 


Feedforward Path 


Hysteresis 

Current Regulator Modulation 

1 




+ ' 

A 





s 


d 


Figure II.4: Power Supply Control 

The nonlinear controller allows for voltage and current regulation as it 
calculates the control signal d that determines the on/off state of the switch. For the 
average value model, d is the duty cycle of the input state. The measured output voltage 
v^ut is compared with the reference voltage and the error signal is fed through a 

proportional plus integral (PI) controller. To provide a fast response to changes in load, 
the load current is added to the output of the PI controller and is multiplied by the 

output of nonlinear stabilizing control (NSC), which has a scaling effect on the influence 
of the output voltage and current. The resulting reference current signal is bounded and 
compared with the measured current in the current regulator. The output of the 

current regulator is used for the hysteresis modulator so that the actual current closely 
tracks the measured current. 
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Parameter 

Description 

Value 

Units 

K./ 

Reference voltage (command output voltage) 

500.0 

V 

i 

max 

Maximum command current 

40.0 

A 

Kr 

PI controller, proportional gain 

1.0 


K, 

PI controller, integral gain 

100.0 


K., 

Current regulator, integral gain 

100.0 


n 

NSC exponent 

1.0 


Krr 

Hysteresis error bandwidth 

1.0 

A 

^Lmax 

PI controller, integral gain output max 

100.0 



PI controller, integral gain output min 

-100.0 


tf 

J max 

NSC Transfer Function output max 

1000.0 


tf 

d min 

NSC Transfer Function output min 

400.0 



Table II.2: Power Supply Control Parameters 


d. Nonlinear State Space Model 

Looking at the power supply average value model circuit in Figure 11.3, 
the inputs to the system can be defined as 


Ml 


4iE 

M2 



M3 


dv^ 

U^\ 


_ _ 


16 


















































If the controller is considered external to the model, the control signal d 
can be viewed as an external input dynamically coupled to states of the system. In this 
manner, nonlinear terms are defined as external inputs and Mj . The model is then 

expressed as a linear state space model driven by inputs Kj and with a decoupled 
nonlinear term as shown 


where 



u, 


U. 


1 

+ ^2,ps 

L 

U,j 



(0.4) 




-[R,^ + 3coJ^l7r) -1 


L. +2L 

dc c 

1 


L, +2L 

dc C 


0 0 

0 0 


0 

0 

-R 




out 

1 




= 


3^/3 

0 


■ 0 

0 ■ 



-1 

0 


0 

0 

, B , = 

~ C , 



0 

0 

’ 2,p.v 

0 

1 

0 

-1 


^out 

0 _ 

_ 0 


c 

out 



and the output equations are chosen as 




1 

0 

n 

1 

0 

0 

0 

0 


0 

1 


0 

0 

0 

0 

0 


0 

0 


0 0 


0 

1 

0 

-R. 


0" 











~0 

o' 



“0 

o' 
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0 

0 



0 
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0 


1 

0 



0 

0 


0 


0 

0 

Kj 


-1 

0 

U 2 

x + 

0 

0 

M4 

+ 

0 

0 


0 




- 




1 


0 

0 



0 

0 


0 


0 

0 



0 

1 



0 

-1 



0 

0 


0 










(0.5) 
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Although the dynamic model has only four states, we included a larger 
number of observations for the purpose of redundancy and fault detection. Figure 11.5 
shows a simulation of the power supply, plotting the input and output voltage for both the 
detailed and average value models. Therein, the input voltage is shown to have the same 
dc value between the two models, but the detailed model shows a small amount of ac 
content and switching dynamics. The output, however, is nearly identical. 


Power Supply Input: 



Power Supply Output: 



Figure II.5: Comparison of Detailed and Average Value Models of 

the Power Supply 


2. Bus 

The bus model connects the power supply to the SSCM in each zone. It is a 
medium length Nominal-T circuit transmission model with switches separating the three 
zones as shown in Figure II.6. The series impedance and shunt admittance parameters 
are equivalent in each zone and are shown in Table II.3. 
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Zone 1 Zone 2 Zone 3 



Figure II.6: Bus Model 


Parameter 

Description 

Value 

Units 

R 

Bus line resistance 

34.0 

mCl 

L 

Bus line inductance 

2.26 

//H 

C 

Bus shunt line capacitance 

0.5 

//F 

G 

Bus shunt line admittance 

10.0 

Q ' 

Rs 

Bus zone short line resistance 

1.0 

Q 


Table II.3: Bus Model Parameters 


Included in the bus model is a shunt resistance used to model bus shorts. 
During no-fault operations, the shunt resistance is infinite to model an open line. While 
not part of the power electronics modules, the bus model has proven to be quite an 
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important link between modules. Generalized forms of the bus model are developed such 
that n zones can be considered on larger scale systems. 


a. Six-state Bus Model 

Consider three operating modes of the bus system where Mode 1 
represents the power supply providing power to Zone 1 only, Mode 2 is where Zone 1 
and Zone 2 are energized, and Mode 3 is where all 3 zones are energized. The mode 
switch configurations for each zone are described in Table II.4. 



Switch 

Switch ^2 

Mode 1 

Open 

Closed 

Mode2 

Closed 

Open 

Mode 3 

Closed 

Closed 


Table II.4: Bus Zone Switch Configuration 


It is desired to express the bus model in matrix equation form. From 
Kirchhoff’s conservation laws, the following equations are derived from the circuit in 
Figure II.6: 


R 1a 1 e 

x = --x + jQ^v^+-Su 


K ~~QLx - v-Mu 

Z ^ — 


= —l+Z 


C " 


( 0 . 6 ) 

(0.7) 

( 0 . 8 ) 

(0.9) 
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where 


_ 


hzy 



U2 

_ 

^bZ2 

M3 


^hzi, 

_M 4 _ 





ps 


x = [x^ X 2 -^3 = bus current through zones , 

z = [zi ^2 ^3 ]^ = voltage across zone shunt capacitance, 

Vj. = J = zone bus voltage, and 

i, = s zone shunt admittance line current. 

Z 1 _ Z] Z2 Z3 J 


From equations (0.6) through (0.9), it can be shown that the six-state 
model can be expressed as 


R 

L 


( 1 ~ 

G Xrr. 


1 

r 1 ~ ^ 

^ QmQm 


p _ 




y RG j 

LG 

X 


L 

y G ” j 

G -J- 

G , 

_z_ 



G 

- —T Q 

- —I 




.M 

CG " ” 

RsCG _ 



CG 


( 0 . 10 ) 


with output equations 
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- 

X 


h 



06x4 

z 

= 

KJ KJ 

1 1 

1 _ 1 

-1- 

G 

h\ 


G Rfi 



\ 

1 

1 _ 


( 0 . 11 ) 


where 


G = 



( 0 . 12 ) 


M = 


1 

0 

0 


0 0 
1 0 
0 1 


0 

0 

0 


(0.13) 


and 


S = 


0 

0 

0 


0 0 
0 0 
0 0 


0 

0 

1 


(0.14) 


The effective conductance G shown in (0.12) includes the shunt resistance term used in 
modeling a bus short. Under no-fault conditions, the shunt resistance line in Figure 11.6 is 
modeled as an open circuit (7?^ = oo). The model also includes the three mode 

configurations with the matrix defined as 
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0 0-1 

1-10 ,m = 1 (mode 1) 
[-1 0 0 
To 1 -f 

= < 0 -1 0 , m = 2 (mode 2) 

[-1 0 0 
To 1 -f 

1-10 ,m = 3(mode3) 
[-1 0 0 

Note: = Ql and = Q]. 


State ordering assignments for the port and starboard busses are also 
considered in the model using the transformations 


"1 0 o 1 To 0 r 

=010 and = 010 

0 0 ij [l 0 0 

to account for opposite-side state numbering. To demonstrate the bus model, a baseline 
scenario is identified for all fault diagnosis experiments. In further chapters, any induced 
system faults will reference the baseline system dynamics for analysis. Table II.5 defines 
the nominal plant configuration. 



Table II.5: DC ZEDS Baseline Scenario Parameters 
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Herein, both busses are in Mode 3 with all six SSCMs on at the start of the 
simulation. Each zone supplies an identical constant power load of 5.65 kW. From the 
Oring Function dynamics, current is supplied from the starboard bus in zone 1 and from 
the port bus in zones 2 and 3. Figures II.7 - II.9 show an example of the bus current 
distribution of the DC ZEDS baseline scenario for both the detailed and average value 
models. 


c 

O) 


3 

o 


Starboard Bus Current. Zone 1 



Port Bus Current, Zone 1 



Figure II.7: DC ZEDS Baseline Scenario Bus Current Distribution, Zone 1 


From the results in Figure II.7, it can be seen that the load in each zone 
requires about 12 Amps. For Zone 1, the port bus (bottom) is directly connected to the 
port power supply, which shows current flowing through to Zone 2 at t=0.3 seconds. An 
additional 12 Amps is added to the port bus, which is being sent to Zone 3 at t=0.4 sec. 
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Current (A) Current (A) 


The starboard bus (top) is at the end of the line since it is directly connected to Zone 3 at 
the opposite end. Here, Figures 11.7 - II.9 show that the starboard bus only feeds Zone 1. 


Starboard Bus Current. Zone 2 



Port Bus Current, Zone 2 



Figure II.8: DC ZEDS Baseline Scenario Bus Current Distribution, Zone 2 
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Starboard Bus Current, Zone 3 



Time (sec) 

Figure II.9: DC ZEDS Baseline Scenario Bus Current Distribution, Zone 3 


b. Three-state Bus Model 

The three-state bus model uses the same parameter values as the six-state 
model as shown in Table 11.3 with the exception of completely passive line impedance in 
each zone. This is equivalent to setting the line inductance L to zero. From a modeling 
and simulation point of view, the states to describe the dynamics of the bus current can 
no longer be found by integration, but by algebraic relationship. The six-state model is 
essentially a double-cascaded, medium length Normal-T circuit transmission line model 
described in [323]. The energy storage elements are very small and equate to transient 
dynamics three orders of magnitude faster than that of the power electronics. It also 
implies that most of the power is transferred from the power supply to the SSCMs with 
very little energy storage in the bus, which is desired. Reducing the number of states may 
improve the realism of the model, but might incur unresolved algebraic loops in the 
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model. It will be shown that a three-state model can be derived and is demonstrated to be 
useful for input estimation on the unknown input observers. 

To begin the model derivation, the same matrix definitions are used except 
this time the three modes are expressed as subspace projections of the full mode. 
Consider the bus model in Figure II.6 and the switching configurations in Table II.4. If 
the bus is operating in the full mode (Mode 3), it can be represented as an orthogonal 
basis in three dimensions as 


/j, 


with each column (x,y,z from left to right) representing the bus current through Zone 3, 
Zone 2, and Zone 1 (from left to right) respectively. The numbering convention in 
reverse order is an unfortunate result of having 2 power supplies entering the DC ZEDS 
system from opposite zones. Each lesser mode can be expressed as a subspace basis of 
<^ 3 . Eor example, the Mode 2 basis spans the y-z plane as 


‘?2 


0 

1 

0 


0 

0 

1 


while the Mode 1 basis spans the x-z plane as 


= 


1 

0 

0 


0 

0 

1 
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In a physical sense, the Mode 2 basis describes opening the switch 
between Zone 2 and Zone 3, which sets the current in Zone 3 to zero and the dynamics 
follow. Likewise the Mode 1 basis is like opening the switch between Zone 1 and Zone 
2, which sets the current in Zone 2 to zero. Now define the mode projection matrix as 


= 



q 


T 

m 


where m = 1, 2, 3 denotes the mode. Removing the inductance in the bus line changes 
the dynamics of the bus and shunt line current Kirchhoff law equations to 


X = 


~Qy, +—Su 
R " R - 


(0.15) 


and 


= -Qx - Mu . 


(0.16) 


With some substitution and algebra on equations (0.8), (0.9), (0.15) and (0.16), the three- 
state model is expressed as 


z =-- 


CG 


T,QC,+j-I 


CG 


{T^QD,+M)u 


(0.17) 


Where G , G , T^, C, R , R^, and M are the same as in the six-state model and 
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Q = 


0 

1 

-1 


1 

-1 

0 


-1 

0 

0 


The corresponding output equations (defined similarly to the six-state model) are 


X 





z 



1+ 

A 



C 3 


A 

J‘z_ 




_^ 4 _ 


( 0 . 18 ) 


where 




RG 


( 1 2 , 
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R 
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A = - 

' R 

D ,=0 


1 


/ + 

V RG 


P 


D,=-Ut,QD,+M) 

Lr 

d, = -^{t,qd,+m). 

u 



\ 


The properties of the 3-state model show that the D matrix is full rank, 
which turns out to be important for full input estimates from unknown input observers as 
will be explained in the following chapter. 

3. Ship’s Service Converter Module 

The purpose of the Ship’s Service Converter Module (SSCM) is to step the bus 
voltage down to an acceptable power quality and rating for zonal loads. As shown in 
Figure II. 1, the DC ZEDS system has 2 SSCMs in each zone - one receiving input power 
from the port bus, and the other from the starboard. All six SSCMs are identical in 
topology and parameters. Rated at 8 kW, the SSCM accepts a 500 V input (that may 
vary) and provides a clean, 400 V output and a maximum of 20 Amps. 


a. Detailed Model 

The SSCM topology is a standard Buck Converter as shown in Figure 
II. 10 with parameters listed in Table II.6. The signal d for the detailed model controls 
the semiconductor switch position to achieve the desired command voltage. 
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Figure II.IO: Ship's Service Converter Module Circuit 


Parameter 

Description 

Value 

Units 

h- 

Input filter inductor 

357.0 

//H 

RfL 

Input inductor resistance 

0.2 

Q 


Input filter capacitor 1 

500.0 

//F 

^fC 

Input filter series resistor 

1.0 

Q 


Input filter capacitor 2 

45.0 

//F 

^Ollt 

Buck Converter output capacitor 

500.0 

//F 

^oul 

Buck Converter dc inductor 

3.0 

mH 

K. 

Resistance of the dc inductor 

0.5 

Q 

Pout 

Rated power output 

8.0 

kW 


Table IL6: SSCM Parameters 
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b. Average Value Model 

Similar to the power supply model, the SSCM average value model 
simulates a fast-average of the semiconductor switching dynamics of the detailed model 
such that state variables are constant in steady state. Again, the average value model 
maintains transient response accuracy and implements the same model and control 
parameters. The model topology depicted in Figure II. 11 shows the averaged command 
signal (duty cycle) d being multiplied by states of the system, which acts as a dependent 
voltage and current source. 


m 



1 


Load 



Figure 11.11: SSCM Average Value Model Topology 

c. SSCM Control 

Because of the fast semiconductor switching action of the SSCM and 
near-ideal nature of the control scheme, the desired output voltage can be met even in the 
presence of input disturbances. For this DC ZEDS model, the controller design is 
depicted in Figure 11.12 with associated parameters listed in Table II.7. 
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Figure 11.12: SSCM Control Architecture 


Parameter 

Description 

Value 

Units 

K./ 

Output reference voltage 

400.0 

V 

K, 

PI controller, proportional gain 

0.1 



PI controller, integral gain 

23.0 


h, 

Current feedforward proportional gain 

0.015 


f 

J sw 

switching frequency 

20.0 

kHz 


Table II.7: SSCM Control Parameters 


Here, the output voltage is regulated through a PI controller design (with 
gains and K.). For simultaneous current regulation the error signal, based on the 

difference between 4 ^, and and scaled by h ., is placed on the feedforward path and 
added to the PI control output. 
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d. Nonlinear State Space Model 

Similar to the power supply model, the control signal d is viewed as 
external to the model and dynamically coupled to states of the system such that the inputs 
can be defined as 


Ml 


1 

s; 

1__ 

M2 



M3 





1 

_1 


In the same manner as the power supply model, nonlinear terms are 
expressed as unknown inputs while the remaining linear dynamics can be viewed as 2 
subsystems (pre-filter and buck portions) coupled together as shown: 


^ ^sscm^ ^l,sscm 


+ B. 


2,sscm 


(0.19) 


with plant dynamics 
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and known and unknown input matrices associated with linear and nonlinear inputs 
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The states of the SSCM depicted in Figures 11.10 and 11.11 are the energy storage 
elements from left to right in the model as 


X = 


c/l 


"c/2 


It should be noted here that there is a limit to the direction of the current 
through the inductor path such that is always in the positive direction defined in 

Figure 11.10 and Figure 11.11. The model handles this state restriction by setting a lower 
bound on the integrator to zero. The output equations are defined as 
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( 0 . 20 ) 


Notice that the feed-through term in (0.20), defined as the D matrix, is 
full rank, which implies that full input estimation can be achieved from an unknown input 
observer. The SSCM model will be used as an example in the following chapter to 
demonstrate the unknown input observer methods with input reconstruction. Simulation 
results from the Zone 1, Starboard SSCM are shown in Figure 11.13 where the DC ZEDS 
baseline scenario depicted in Table II.5 is used. 
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SSCM Stbd Zone 1 Input: V. 



Figure 11.13: Zone 1 Starboard SSCM Simulation Regulates 400 V Output (bottom) 

from 500 V Input (top) 


4. Oring Function 

The Oring function is the model in Figure 11.1 that determines which bus supplies 
the power to the load in each zone. The port and starboard SSCMs have slightly different 
voltage regulation set points so that the SSCM output voltage with the higher set point 
picks up all of the required power to the load in the zone. The set points are close so that 
if the primary SSCM is lost, then the alternate SSCM picks up the load requirement at 
nearly the correct voltage. As an example, for the nominal simulation configuration 
shown in Table II. 5, Zone I load current is supplied from the starboard bus as shown in 
Figure 11.14. 
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Load.s 



Time (sec) 

Figure 11.14: Oring Function Simulation Results 


The function is modeled as two voltage sources on opposite ends of a circuit 
providing current to a common constant power load through opposing diodes as shown in 
Figure 11.15 and parameters detailed in Table 11.8. If the inputs and states of the state 
space model are defined as (respectively) 



r ' 





u = 

V, 

and 



^Load,p 


^load _ 




^Load,s 


then the state space model is defined as 
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The opposing diodes in the model are simulated as lower bounds on the state integrators. 


1 


Load.p 


V 

^ Load 


^ Load. 



Figure 11.15: Oring Function Circuit Model 
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Parameter 

Description 

Value 

Units 

Kr 

Port/stbd line resistance 

1.0 

mCl 


Port/stbd line inductance 

100.0 

//H 


Table II.8: Oring Function Parameters 


5. Inverter Module and Constant Power Load 

According to the DC ZEDS model, the load in each zone is a Ship’s Service 
Inverter Module - converting dc voltage back to ac - followed by a constant power load 
bank. Although the simulation model is more detailed, the relevant dynamics may be 
represented by a capacitor and effective series resistance in parallel with an ideal constant 
power load of 5.65 kW. Herein, the EDI analysis is conducted on the DC side of the 
power distribution system, which encompasses the rectifier dc input voltage of the Power 
Supply to the dc voltage and current leading to the SSIM in each zone. 


C. METHOD FOR ANALYZING DC ZEDS SYSTEM 

With respect to the complete set of state space models, including power supply 
and SSCM control, the DC ZEDS model is comprised of 21 distinct state space models 
with a total of 83 inputs, 159 outputs and 68 states. Among the inputs, 2 are external plus 
8 reference voltage set points. The remaining 88% of inputs are either output 
measurements from connecting modules or nonlinear components that are a function of 
control signals and internal states. This fact demonstrates the high degree of dynamic 
interdependence between the modules and also the model’s propensity to propagate faults 
from one module to another. Among all the output measurements, only about 30% of 
them are directly used as inputs to controllers or other DC ZEDS modules. The 
positioning of the remaining 70% of output measurements is based solely on enabling 
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observer sufficiency conditions for full-state and input estimation. In other words, almost 
70% of the sensors throughout the DC ZEDS system are used exclusively for detecting 
faults and analytical redundancy. 

It is shown that the average value model accurately represents the signature of the 
detailed model - even in the transient response. In steady state, a frequency analysis 
between the two models reveals additional frequency content in the detailed model due to 
the switching action and it can be made to closely resemble the AVM dynamics with a 
simple lowpass filter (cutoff frequency of 100 Hz). 
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III. A POLYNOMIAL APPROACH TO DYNAMIC MODELING 

AND FAULT DETECTION 


The fundamental framework to detecting fault characteristics in this work is based 
on analytical redundancy. Indeed, the large-scale dynamically interdependent DC ZEDS 
system must rely on mathematical models to provide fault indications for two main 
reasons: first, as outlined in Chapter I, it is efficient in terms of hardware invested 
compared to current methods [390][271][272][273][289][343][90]. The second reason is 
that the DC ZEDS model is well-defined as shown in Chapter II, which lends itself well 
to model-based estimation techniques[220][335]. By combining the physics-based 
dynamical models and sensor information, we compute estimates of the states of the DC 
ZEDS modules. Then, a comparison of the sensor measurements and their estimates is 
made, yielding residuals that are used as a measure of input-output consistency. 

Perhaps the most widely used approach to generate residuals is to measure the 
consistency between observations and their prediction based on a mathematical model 
[60]. Standard approaches are the observer for deterministic linear and nonlinear models 
and the Kalman Eilter for stochastic models. In both cases, the model-based estimation 
technique require knowledge of all inputs u{t) and outputs y{t). There are cases 
however (especially in residual generation), where it is desirable to estimate the system 
when one or more (perhaps even all) inputs are unknown. Eor example, it could be that 
the dynamical model of a system can be separated into linear and nonlinear parts and the 
nonlinear part could be modeled as an unknown input. Also, it might be desirable to 
remove the dynamical influence of a particular input from an observer estimate for the 
purpose of detecting and isolating actuator faults. In either case, it has been shown that 
unknown input observers are well-suited to solve the problem [110] [175] [297]. 

One such method to obtain an unknown input observer is the subspace method 
where the full state is estimated based on dynamics residing in a subspace of known 
inputs and outputs as stated in [170] and [175]. The conditions for UIO existence are 
proved in [60] and [230] in which a full-state unknown input observer is assumed. 

To this goal, consider a dynamic system 
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J x{t) = Ox + FjMj (0 + T 2 U 2 (t) 
\y{t) = Cx{t) 


and we want to find conditions to estimate the state x without the input u^it). This 
problem can be solved by decomposing the state in two components [175] 

X = (/ -MC)x + MCx 
= q + My. 


The matrix M is chosen so that -MC) is a projection matrix and q = (/-MC)x is 
independent on the input u^it). The component q(t) satisfies the differential equation 


q{t) = {l-MC)^x{t) + {l-MC)T^u,{t) + {l-MC)T2U2{t) (0.23) 


where Mj( 0 and u^it) represents all the known and unknown inputs, respectively. 
If M is such that 


{1-MC)T2=0 


(0.24) 


q{t) is clearly independent of u^it). For a more detailed explanation of the subspace 
method, refer to Appendix C. 


Lemma 4.1 Equation (0.24) is solvable iff: 
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mnk{CT ^) = rank (Fj) 


(0.25) 


Proof: Necessity : When (0.24) has a solution M , it becomes MCT^ = Fj or 

(cr,f H^=r/ 

which means F 2 belongs to the range space of the matrix (CF 2 )^ resulting in: 

rank(T2^< rank^i^CT^Y^ 
and 

rank (F2) < rank (CF2). 

However, 

ranA:(CF2) < min|ranA:(C),ranA:(F2)} < rank(T ) 

Hence, rank (F 2 ) = rank (CF 2 ) and the necessary condition is proved. 

Sufficiency : When ranki^CT^) = rank(T 2 ) holds true, CF 2 is a full column rank 
matrix (because Fj is assumed to be full column rank), and a left inverse of CF 2 exists: 

(cr,)’=((cr,f 

Clearly, M = T 2 (^ 2 )^ is a solution to (0.24) [60]. The rank condition also implies that 

the number of independent output measurements must be greater than the number of 
unknown inputs [175]. 

A full-state UIO estimate is achievable if the rank condition (0.25) is satisfied. 
This condition can be rather restrictive. A simple example of a system in controllable 
canonical form would show that a necessary condition would be that at least one of the 
transfer functions to have a relative degree at most equal to one. 
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For the purpose of fault diagnosis, what is actually needed is a residual generator, 
not necessarily the full-state estimate. In the case where the full-state estimate is not 
achievable, it is still possible to compute the residual based on the partial state. The 
purpose of the polynomial method is to take advantage of this property by generating 
residuals based on the partial state. 

The remainder of the chapter develops two unknown input observer (UIO) 
methods generalized to the multivariable case on a polynomial approach. The first 
process, called the input replacement method, is limited to a more restrictive class of 
systems, but has a simpler combined input and state Kalman estimation implementation 
algorithm suitable for fault diagnosis. By comparison, the second method (called the 
parametrization method), takes a more general class of systems and is based on the 
nonuniqueness property of the Diophantine equation to establish an unknown input 
Kalman estimation implementation. In both cases, a partial-state representation of the 
system is sought, which is directly related to the controllable canonical form. 

A. DYNAMIC SYSTEMS AND THE PARTIAL STATE 

I. Single Input Single Output (SISO) Systems 

Consider a single input single output (SISO) system represented as a proper 
rational fraction of scalar polynomials B{p) and A{p) as shown in Figure III.l and 
expressed as a transfer function along with its relationship to the state space realization in 
( 0 . 26 ). 


Hit) 

B{p) 

y(0 


Mp) 



Figure III.l: SISO System in Transfer Function Form 
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H{p)=^^ = C{pl-(^X'T + D 
A{p) ’ 


(0.26) 


The numerator is of order or at least of order -1 = O) for a proper 

or strictly proper transfer function/f(p), respectively, and is of the form 

B(p) = b^p"' + b^p ' + b^p'''^^ +... + b^^. (0.27) 

Likewise, let the denominator characteristic equation be a monic polynomial of 
order as 


A(p) = p"'+a,p"'-‘+... + «„ (0.28) 

where p is the differential operator s in continuous time and the time-shift operator z 
in discrete time. If there are no common factors (i.e. no pole-zero cancellations) between 
B{p) and A{p ), then the transfer function is said to be coprime and the system is a 
minimal realization. Now referring to Figure III.2, define z{t) as the partial state [200] 
(or sometimes called the pseudo state [61]). 



Figure III.2: Partial-state Space Representation of SISO System 
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It is easy to see that it can be related to the input and output signal, u{t) and y{t) 


as 


A(p)z(t) = u(t) (0.29) 

y(t) = B(p)z(t) (0.30) 


The meaning of the partial state is that (0.26) can be related directly to the state of a 
controllable canonical form realization. Again, in the SISO case this is straightforward, 
since the state of the system is 


z(t) = [z(t), pz(t),..., p"‘ ‘z(0] 


(0.31) 


where z(t) is a scalar and z(t) e M"’. 


2. Multiple Input Multiple Output (MIMO) Systems 

Expanding the discussion of partial state to the more general multivariable case, 
consider the class of systems represented in state space as 


px(t) = (l>x(0 + Tu(t) (0.32) 

yit) = Cx(t) + Du(t) (0.33) 


where veM"' , mgM"", y and 0,r,C, andD are of appropriate dimensions. 

Again p in the state dynamics equation (0.32) denotes the differential operator in 
continuous time or the time-shift operator in discrete time and its notation is consistent 
throughout. Equations (0.32)-(0.33) can also be represented in transfer function form. 
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only the “numerator” B{p) and “denominator” A{p) are matrices whose elements are 
polynomials. Together, H{p) is called the transfer function matrix or simply the 
transfer matrix and can be written as 


H{p)^B{p)A-\p) (0.34) 

H{p)^A-\p)B{p) (0.35) 


The right fraction pair B{p) and A{p) has dimensions n^xn^ and x , respectively 
while the left fraction pair A(p) and B(p) has dimensions n^xn^ and n^xn^. The 
“denominator” polynomial matrices A{p) and A{p) are always square. If the pair 
(0,r) is controllable and (©jC) observable then the system is a minima! realization, 

and the polynomial pairs (5(p),A(p)) and {^B(p),A(p)^ are mutually coprime. In the 

SISO case, this implies that the numerator and denominator polynomials have no 
common roots (no pole-zero cancellations). The more general MIMO case requires an 
appropriate extension of this concept, as shown in Appendix D. Under these conditions, 
the three representations 


((b,r,c,D) 


in state space form. 


{A(p),B(p)) 


in right polynomial matrix fraction form, and 
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{A(p),B(p)) 


in left polynomial matrix fraction form are equivalent and one can be computed from the 
other [61]. Of particular interest is the right coprime factorization in (0.34) which leads 
to the MIMO partial-state representation of the system shown in Figure III.3. 



Figure III.3: Right Fraction Representation of MIMO System 

By extension to the SISO systems, the partial state in the multivariable case is related to 
the input and output signals u(t) and y(t) of (0.29) and (0.30) and the system can be 
realized as 


p'^'Ziit) 


-VAi(0 + 4^‘w(0 


y{t) = (S, - B.A;'A)i«) + BXm 


(0.36) 


(0.37) 


where the state space equations (0.36)-(0.37) are related to the original system (0.32)- 
(0.33) by a nonsingular matrix T such that 


x{t) = Tz{t). 


(0.38) 


50 




The matrices are matrices of high degree coefficients whose dimensions are 

n^xn^ and n^,xn^, respectively, and A,, 5, are the matrices of all lower order 
coefficients with dimension x and n^xn^, respectively. It should also be noted here 

that a method for determining T in (0.38) to transform the original system in (0.32)- 
(0.33) directly to control canonical form in (0.36)-(0.37) is found in [386]. A more 
detailed explanation of the system partial-state realization for the MIMO case is found in 
[61] and in Appendix D. 

Similar to the vector of partial states (with derivatives) in (0.31) for the SISO 
case, we see that the partial state is a vector related to the system as 


z(t) = 


z, (t). 


* 7. ifi • 


(0.39) 


where z(0el^”' and called the column degrees (i.e. the 

maximum degree of the polynomials in each column) of the denominator matrix A(p) 
and is related to the controllability index in [200] and [386]. From (0.39), it can be seen 
that the maximum degree of the polynomials in each column of the denominator 
polynomial matrix A{p) is directly related to the full-state estimate of the observer. It is 
an important distinction, then, to define the partial state alone (without derivatives) as 


z(0 = [zi(0,Z2(0,"-,z„_^(0] 


(0.40) 


where z{t) e M"" has exactly the same elements in z(t), only without the associated 
derivatives. This reduced partial-state vector is important as it relates to developing 
unknown input observers when the UIO rank condition is not satisfied. 
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It should be noted here that in Fault Diagnosis, the goal is to generate residuals, 
not necessarily to estimate the full state. We will see that in many situations, it is 
possible to estimate the partial state z{t) but not the entire state x{t) = Tz{t ), especially 
when u(t) is unknown or assumed unknown. In this case, however, if a full-state 
estimate must be made, it still can by including the appropriate number of derivatives to 
the partial-state vector through a variety of filtering methods to obtain z(t) in (0.39). 


B. POLYNOMIAL-BASED UIO BY INPUT REPLACEMENT 
I. Theory 

In order to illustrate the UIO concept, let us consider a simple case of a SIMO 
system with equations (0.32)-(0.33) with one input and two outputs as shown in Figure 
III.4. The goal is to estimate the partial state z(t) based on outputs y^(t) and y 2 (t) only. 



Blip) 

}\i') 

1,(0 


Aip) 






B.ip) 

yA’) 


Aip) 



Figure III.4: SIMO System Transfer Function 


As previously discussed, z(t) satisfies the equations 

Mp)z(t) = u(t) 
y lit) = Blip) z(t) 
yiit) = B2ip)zit). 
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If the input u{t) is not available, we can use the last two equations as 


B^{p)z{t) = y^{t) 

yiit) = B^{p)z{t) 


with the same partial state z{t). 

From the definition of the partial-state vector in (0.31), a set of equations can be 
derived from the input/output relationships as if y^it) and jjCO the input and output 
of the system as shown in Figure III. 5. 



Figure III.5: SIMO Unknown Input System 

Figure III.5 shows that z{t) is the partial state associated to the transfer function 
B^{p) jB^{p) and therefore it can be estimated by a standard observer or Kalman Filter. 

In order for the UIO to exist, B^{p) and B^{p) must be mutually coprime (no 

pole-zero cancellations), which is the same as saying the system in Figure III.5 must be 
observable. Note that there are more linearly independent outputs (2) than there are 
unknown inputs (1), consistent with the rank condition (0.25) [230]. To add robustness 
to the state estimation, a standard Kalman Filter can be used to minimize the mean square 
error of the estimates. 
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In general, consider a MIMO system of the form (0.32)-(0.33) where at least one 
of inputs is unknown and the input matrix T with corresponding input vector function 
u are partitioned into known and unknown inputs as 


rM(0 = [r, 



(0.41) 


where input is partitioned into known inputs and =n^ unknown inputs. In 
this manner, using equation (0.41) we can write 


x(t) = Ox(0 + FjMj (0 + FjMj (0 (0.42) 

y(t) = Cx(t) + Duit) (0.43) 

where Fj e and Fj e . For the unknown input observer to exist for full-state 
estimation, rank (CF 2 ) = rank (F 2 ), which implies that there must be at least as many 

independent outputs as there are unknown inputs. From the right polynomial matrix 
fraction of the transfer function matrix H(p), combine the partial-state input/output 
equations (0.29) and (0.30) as 


A(p) 


u(t) 


z{t) = 


B(p)j 




(0.44) 


where A(p) e , B(p) e , z(t) e R"", u(t) g R"“ , and y(t) g R"^ . It can be 
seen in equation (0.44) that there are a total of +n^, equations relating the input and 
output measurements to the partial state z(t ). The main idea is to replace all the 
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unknown input row equations in A{p) with the same number of output row equations in 
B{p) to obtain a new transfer matrix 


H(p) = B(p)A-\p) 


(0.45) 


where "" " and A(p) e . Right coprimeness of H(p) and column 

reducedness of A(p) must be verified since there is no guarantee that the system will 
have satisfied those conditions. 

To illustrate this argument, suppose a system described in (0.42)-(0.43) where 


= 5, = 2, and n^=4 and suppose 

equation (0.44) is represented as 


AoP^ + + Ai : 

l^loP 2^\lP l^llP 2^13 


AoP^ + + AiP' + 1^13 

l^wP^ + l^nP^ + l^ilP + 2^13 
3^ioP AlP i^nP An 
AwP^ + Aiy + AnP + 4^13 


the second input u^it) is unknown. Then 


I^IqP l^llP 1^22 


uA) 

2 1 

2^2oP 2^2\P 2^22 


U^it) 



zA) 



1^20P" + 1^2lP‘ + iK 




yA) 

2^2oP 2^2\P 2^22 


1 

to 

1_ 


3 ^ 2(0 

3^2oP ■*" 3^2lP ■*" 3^22 


3 ^ 3(0 

4^20P" + 4^221^* + A22 _ 




It is desired to have the third output y^) replace the unknown input u^it) in the input 
polynomial matrix A(p). The rank condition has been satisfied and assuming the 
leading coefficients on the polynomials for the third output, and jZjjq are nonzero 
(which implies that in (0.36)-(0.37) is nonsingular), the replacement can be made as 
shown 
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\^\oP ^ l^llP l^llP 1^13 

i^loP 3^l\P i^llP 3^13 • 


AoP'+ Ay + AlP + A 3 '■ 
iKP^ + l^y + ihlP' + 2^13 '■ 

aKp^ + Ay + AlP + A3 '■ 


l^2oP \^2\P 1^22 


1 

1_— 

3^2oP 3^2lP 3^22 


>^ 3(0 



zA) 



Ay + 1^2iF' + A 2 



= 

yA) 

2^2oP 2^2\P 2^22 




3 ^ 2(0 

A^y + aKA + A22_ 


1 

1 _ 


The top = 2 rows are now the new input polynomial matrix A(p) and the bottom 
-1 = 3 rows comprise the new B(p). The transfer matrix H(p) = B(p)A \p) is 

checked for right coprimeness and A(p) is verified to be column reduced. Finally, the 
new system of equations is represented in state space form for Kalman Filtering 
implementation. Note that the estimate for the new system will be on the partial state, so 
it must be transformed through the linear transformation in (0.38) if the original full-state 
estimation is desired. 

The relationship of the UIO rank condition in (0.25) to the polynomial 
representation for full-state estimation implies that B(p) must have higher degree 
coefficients. Consider the input and output matrices F,C for the controllable canonical 
form realization both partitioned as 

- d 

c =[c, - c._ 


(0.46) 


with 


Ff =[0---0r,.]Gi?"""« 


(0.47) 
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for i = It turns out that the terms C; are the coefficients of the polynomial 

B(p). In particular we can write 


A(p) = 

B(p) = 


n 


Yn. 


E (p) + terms of lower degree 


(0.48) 




CiP^ ••• 


-1 


+ terms of lower degree 


where 


(p) = diag^p'^',p^\..., p^"" I e 


In the case of a UIO we can see some of the conditions of existence. Let L be 
partitioned as L = [rj Lj] where r 2 is associated with the unknown inputs. In Kudva 
et al [230], the condition for existence of the UIO for full-state estimation is 
rank(CLj) = rank(r 2 ). From (0.46)-(0.48) it can be seen that 

(0.49) 

(=1 

with the part of the vector associated to the unknown input. From (0.48) we see 
that a necessary condition is that not all p can be zero, which means that the polynomial 
B{p) must have higher degree terms. 

While the unknown input observer rank condition (0.25) must be satisfied for full- 
state estimation, it does not need to be satisfied for estimating the partial state. This is the 
advantage of using the polynomial UIO technique and is illustrated in the following 
example. 
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2. Example 

Consider a third-order, single-input, two-output system shown in Figure TTT.4 
where the continuous time system is represented in right matrix fraction form as in (0.34) 

with 5 ( 5 ) = [ 5 - 1-1 5-1-2]^ and A(5) = 5 ^- 1 - 65 ^- 1 -II 5 - 1-6 . It can be seen that the 
greatest common right divisor between B(s) and A(5) is 1, which is unimodular. In [61] 
this implies that the system transfer matrix H(p) is controllable and observable and an 
equivalent minimal realization representation of the system in state space form is 



■ 0 

1 

0" 


~ 0 ~ 

x(t) = 

0 

0 

1 

x(t) + 

0 


-6 

-11 

-6 


1 


u(t) 


(0.50) 


y(t) = 


1 1 0 

2 1 0 


x(t) -I- 


u(t). 


(0.51) 


Notice all the terms in (0.50)-(0.51) are directly related to the coefficients of the 
polynomials A(s), and ^ 2 ( 5 ). By observation, the system of equations is indeed 

controllable because it is realized in controllable canonical form where the state vector 
x(t) happens to be the partial-state vector z (0 with derivatives up to order 2 as shown in 
(0.31). Since there is only one input, there is only one partial state z(t). It can be 
verified that the system of equations is also observable. To design a full-state unknown 
input observer, the rank condition (0.25) must be satisfied. For this example, it means 
that 



r 



"o' 

3 


r 

“O' 




“1 

1 o' 







rank 


2 

1 0 

0 


^ rank 


0 



1 



1 



1 

1 

y 
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and a full-state UIO does not exist for this system. However, an unknown input 
observer still exists from a polynomial approach. 

With the polynomial input replacement UIO method described in this chapter, an 
estimate of the partial state can still be achieved as outlined in the following steps to the 
algorithm. We start with the example state space representation in (0.50)-(0.51). 


1) Represent the system of equations in right matrix fraction from as in (0.34) and 
ensure that H(p) is right coprime, which implies a minimum realization. 


H(s) = B(s)A-fs) = 


s + l 
s + 2 


1 


s ^ + 6 s ^ -1-II5-1-6 


The polynomial transfer matrix is right coprime. 


2) Represent the system of input/output equations in partial-state form as in (0.44). 


s^ -\-6s^ -1-II5-1-6 


1 

s 

1 ___ 

5 + 1 

z(t) = 

3^1(0 

5 + 2 


rM 


3) Determine the right matrix fraction form of the UIO where H{p) = B{p)A '^{p) 
by replacing the rows in the input polynomial matrix A(p) that correspond to 
u^lt) with the rows in the output polynomial matrix B(p) that correspond to the 
desired input replacement equations. 
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For this example, there is only one input and two options for replacement. Let us 
choose the output y^{t) to replace u{t). 



5 + 1 


Jxit) 



z(t) = 



5 + 2 


I 

NJ 

1_ 


This leads to the system shown in Figure IIL5 with the new scalar transfer 
function H{p) written as 


H{p) = 


ks) 


B,{s) 

B,{s) 


5 + 2 

7 + 1 ' 


Since H{p) is a first order system - while the original H{p) was third order - 
we can estimate z{t), but not z(t) and z(0 . This means that we can still create a 
residual, but not a full-state estimate. The principle can be generalized to say that 
for any column degree reduction of the original A(5) to the new A(s), it will 
result in an unknown input observer that cannot obtain a full-state estimate. 


4) 


Check H(p) = B{p)A \p) for right coprimeness and column reducedness. 


The new transfer function is 

H{s)^ks)IMs) = '-kk 
5 + 1 
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since B(s) and A(s) are scalar polynomials. By inspection, there are no pole- 
zero cancellations of H(s) and a scalar polynomial is always column (and row) 
reduced [61]. Note that an equivalent observer (in the steady state) can be 
determined from H \s). 


5) Represent H(p) in state space form and implement any standard observer 
method. If the UIO can estimate the partial state (i.e. the UIO rank condition is 
satisfied), determine the nonsingular matrix T such that x(t) = TzXt) to determine 
the original state estimate. 


For this example, the UIO rank condition is not satisfied and a full-state estimate 
cannot be achieved. The state space representation for H(s) is 


z(t) = -z + yi(t) + Vi 

y 2 (t) = z + yft) + v^ 


where the “input” is yft) and the “output” is jjCO • The additional terms Vj, Vj 
account for possible measurement noise. For this problem, let the input 
u(t)= sinf) and the measurement noise variances and be 1% of the 

amplitude of the input signal. Also, all measurement noise are assumed to be 
uncorrelated. Under these conditions, the Kalman Gain is L = 0.5811 and the 
observer implementation equation is 


1(0 =-1.581 lz(0 +[0.4189 0.5811] 


Ji(0 
3^2 (0 


and the output is the partial-state estimate z(t) itself. From the form of the 
observer equation, it can be easily seen that the estimate is independent of u(t). 
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A simulation of the polynomial UIO input replacement method is shown 
in Figure III. 6 and the noise terms were removed in the simulation to see the 
convergence. 



Figure III.6: UIO Input Replacement Method Simulation 


In this example the replacement of the unknown input u{t) with one of the 
outputs y^{t) is straightforward, since the resulting transfer function is 

scalar. In the more general MIMO case, the system resulting from input substitution has 
to be realizable, which might not always be the case. In the next section we present a 
general observer parametrization which has the unknown input observer as a particular 
representation. 
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C. POLYNOMIAL-BASED UIO BY PARAMETRIZATION 

1. Theory 

In the previous section, the UIO input replacement method took advantage of the 
partial-state realization in right polynomial transfer matrix form (0.34). This polynomial 
method, on the other hand, is based on the left polynomial transfer matrix in (0.35). 

In this section a polynomial transfer matrix formulation of the observer or Kalman 
Filter yields a class of stable solutions easily parametrized. Consider the LTI system in 
(0.42)-(0.43) rewritten here as 


x{t) = Ox(0 -I- FjMj (0 + r2U2 (0 (0.52) 

y(t) = Cx(t) + + D 2 U 2 (t) (0.53) 


where xeM"' , mgM"", y and 0 ,r, andC are of appropriate dimensions. 

Furthermore, Fj e M""‘ and F 2 e M”"" are the subdivided constant input matrices such that 
M,(0 and 1 / 2(0 are vectors containing all the known and unknown inputs, respectively. 

If the system (0.52)-(0.53) is a minimal realization, it can be equally realized in partial- 
state form from (0.44) with the additional known and unknown input partitioned as 


u(t) = 


Ml (0 


A(p) 

1/2(0 


_A(p)_ 


z{t) 


y{t) = B(p)z(t) 


(0.54) 

(0.55) 


Now from the general form for the Kalman estimator, a polynomial representation of the 
observer that estimates the partial state z(t) is expressed as 


63 



C{p)z(t) = E(p)u(t) + F(p)y(t) 


( 0 . 56 ) 


with C(p) e , E(p) e M""*'"", and E(p) e . These polynomial matrices can be 
computed either from the Diophantine Equation (Appendix D) 

C(p) = E(p)A(p) + E(p)B(p) (0.57) 

or directly from (0.32)-(0.33) as an observer of the form 

pz(t) = r ' (O - LC)Tz(t) + T-^ru(t) + T-^Ly(t) 

/(t) = Cz(t) 


where y*(t) is the partial-state vector estimate (with derivatives) z(t)eM."‘ in (0.39) 
when the UIO rank condition is satisfied {c = l^ and the partial state z(t) e M"" in 
(0.40) when it is not. In the latter case, C is an x matrix that contains the rows of 
1^ that extracts the partial-state estimate zit) e from the partial-state vector (with 

derivatives) zit) ^ 1^"' • Also, T is the nonsingular matrix transformation to the partial- 
state space shown as x{t) = Tz{t). 

For convergence of the observer in polynomial form, it is shown that the 
substitution of u(t) and y(t) in (0.54)-(0.55) into the polynomial representation of the 
partial-state observer in (0.56), along with the Diophantine equation identity in (0.57), 
yields the state error equation 


c{p){m-z{t))=o. 


(0.59) 
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Equation (0.59) implies that limz(t) - z(t) = 0 provided det(C(p)) has all roots in the 
stable region. 

The useful feature of the Diophantine equation approach is that the solution is not 
unique. In fact, if E^^{p) and Eq(p) e together are a solution of the 

Diophantine Equation (0.57) it can be shown in [61], [200], and [229] that, for any 
arbitrary polynomial matrix Wip) of appropriate dimensions 


Eip) = E,ip)-Wip)Bip) 

(0.60) 

Eip) = E,ip) + Wip)Aip) 

(0.61) 


are solutions to the differential equations. As a consequence, from (0.56) and (0.60)- 
(0.61) we obtain a family of observers, parametrized by iy(p) e as 


Cip)zit) = [E^ip)-Wip)Bip))uit) + [F^ip) + Wip)Aip))yit) (0.62) 


where A(p) e and Bip) e are the left coprime representation of the original 
system transfer matrix Hip) = A \p)Bip). 

Now partitioning uit) into known and unknown inputs we obtain 


Eip)uit) = 


[EM EM]-Wp) ^liP)] 


^niP) ^ 12 ( 7 ’) 

^21 (p) ^22 (p) 


Ml(0 

u^it) 


(0.63) 
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where E,^(p) E’qjCp) e l¥i(p) e W^(p) ^ 

B^,{p) e , 5 i2(p) e , 52i(p) e , and 522(p) e . In order 

to make the observer insensitive to u^it) , we choose the polynomial matrix W{p) so 
that 


E,,(p)-W,(p)B,,(p) + W,(p)B,,(p) = 0. 


It can be seen that the row dimension of B^^ip) requires that , which means that 

in order to have a solution, the number of independent outputs must be greater than the 
number of unknown inputs. This is the same result in Lemma 4.1 where 
mnk{CT^^ = mnk{Tj^. 

If we choose two polynomial matrices W^ip), and W^ip) that satisfy the 
Diophantine equation 


E,,(p) = W,(p)B,,(p) + W,(p)B,,(p) (0.64) 


then a viable solution can be obtained, provided that B^ 2 (p), andB 22 (p) are right 
coprime. Finally, combine (0.62)-(0.64) to obtain 


C(p)z(t) = 


f 

EM-W\p) 

V 


Snip) 

B22iP)_ 


Mi(0+ [FSp)+^\p)^ip)) y(t) 


(0.65) 


with W*(p) = [Wi(p) 1^2(p)] being the combined solution to (0.64). From (0.65), it can 
be seen that the parametrized observer is now independent of 1 / 2(0 • 
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2. Example 

In summary, the steps to obtain a UIO using parametrization is shown in the 
following example. The goal is to demonstrate how the observer is obtained and 
compare the simulation results to the subspace projection method in [175] and to a 
standard Kalman Filter estimation. 

We start with the state space representation in (0.52)-(0.53) where the system, 
introduced in [175], is a fifth-order, two-input, four-output model where 


0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-0.1540 

-0.0042 

1.5400 

0.0000 

0.0000 

0.2490 

-1.0000 

-5.2000 

0.0000 

0.0386 

-0.9960 

-0.0003 

-0.1170 

0.0000 

0.0000 

0.5000 

0.0000 

0.0000 

-0.5000 


■ 0.0000 

0.0000” 










”0 

1 

0 

0 

-1 

-0.7440 

-0.0320 


0 

0 

1 

0 

0 

0.3370 

-1.1200 

, and C = 

0 

0 

0 

1 

0 

0.0200 

0.0000 


1 

0 

0 

0 

0 

0.0000 

O.OOOOJ 








For this problem, both inputs are considered unknown where 
u^it) = [cos(0 sin(t)]^ . Added white Gaussian noise is also imposed on the inputs and 

output measurements where the variance for each is 1.0% of the amplitude of the input. 
All measurements are assumed to be uncorrelated. Given that the initial conditions in 
[175] are arbitrarily chosen to be 


v(0) = [0.3420 0.3200 0.0718 -0.2870 -0.9497]^ 
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and the initial conditions of the observers are set to zero, we set out to determine the 
polynomial-based UIO from the following steps: 


1) Verify the system is a minimal realization and transform the state space system to 
controllable canonical form. 


The system is controllable and observable (and therefore a minimal realization 
[61]). From the system of equations 


zit) = T-^Q>Tzit) + T-^T^u^it) 


y = CTz(t), 


and the system matrices O, Fj, and C previously defined, the nonsingular 
matrix T is determined to be 


-0.0001 

0.0119 

-0.0003 

0.4478 

0.8941 

0.0173 

-0.0009 

0.0000 

-0.0323 

-0.4470 

-0.0012 

-0.1912 

0.0173 

0.0023 

0.2235 

-0.9950 

0.0024 

-0.0001 

0.0894 

-0.0449 

0.0034 

-0.1119 

-0.9950 

-0.0065 

0.0225 


which transforms the original system into controllable canonical form where the 
partial-state vector (with derivatives) is zit) = T 'v(0. 


2) Design a standard observer or Kalman Filter based on the partial-state 
realization. 
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From the model, the plant and measurement noise covariance matrices are set to 
be Q„ = 0.14 K =0.01/„ , respectively. The steady-state Kalman Gain is 

determined to be 


0.0225 

-0.0471 

0.0174 

-0.0135 

-0.0115 

-0.0377 

0.0277 

0.0098 

-0.0513 

0.0563 

-0.0015 

0.0267 

0.0937 

-0.5402 

0.0793 

-0.8628 

0.1225 

-1.3764 

0.3712 

-0.4888 


and the full-state Kalman Filter can be implemented as in equation (0.58). 


3) Check UlO rank condition. Determine C such that 

a. If the rank condition is satisfied, then a full-state estimate can be achieved 
and C = f . 

b. If the rank condition is not satisfied, then determine C such that n^ 
partial states are extracted (no derivatives) from z(t) resulting in the 
partial state z(t) in (0.40). 


The rank condition is satisfied and therefore for this example C = f. 


4) Determine the polynomial matrices C(p), E^fp),and Ff^(p)from (0.51), which 

is the equivalent polynomial representation of the partial-state observer designed 
Step 3. 


In the example, 
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-s-0.0420 0.9300 


C(5) = 


-0.0740 

-s-0.1500 

-1.0000 

-2.2000 

-0.0020 

-0.2500 

12.0000 

26.0000 


0.0720 -0.0150 

0.9600 0.0073 

-s-1.1000 -0.0097 

-0.7900 s-i-0.8700 
3.8000 -0.0560 


-0.0460 

-0.0390 

0.0560 

-0.4700 

-s-2.3000 


and Eq(s) and Fq(s) are constant polynomial matrices 


■ 0.0000 

0.0000" 


■-0.0091 

0.0470 

-0.0170 

0.0130" 

0.0000 

0.0000 


0.0120 

0.0380 

-0.0280 

-0.0098 

-0.1500 

0.0130 

, F,{S) = 

0.0510 

-0.0560 

0.0015 

-0.0270 

0.0000 

0.0000 


0.0940 

-0.5400 

0.0790 

-0.8600 

0.2500 

-I.IOOOJ 


-0.1200 

1.4000 

-0.3700 

0.4900 


5) In (0.62), determine A(p) and B(p) from the left coprime realization of the 
original system, H(p) = A fp)B(p). 


Because the example is a minimal realization, the original system of equations 
defined by the matrices (D , Fj, and C can be realized in left matrix fraction form 
where 


A(5) = 


/+0.8455+ 0.470 
-0.200s-0.098 
0.08S+0.040 
-0.025s-0.012 


2.400s -2.5005-0.610 
-0.0535 + 0.035 
0.5705 + 1.300 
-0.2105 + 0.440 


0.4505 -12.0005-5.100 
0.0175 + 0.029 
0.1305 + 2.900 
-0.0495-1.100 


0 . 955 "-0.415 + 0.010 
-0.0895-0.001 
-0.6805-0.005 
-0.655 + 0.002 


and 
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0.6705-0.280 2.7005 + 1.100 
0.130 0.066 

0.140 -0.640 

-0.054 0.240 

Notice that the sum of the row degrees (maximum polynomial degree of each 
row) of A( 5 ) is equal to the total number of states =5), which implies that 

A( 5 ) is row reduced. Since the degree of the i"" row in B{s) is less than the 
degree of the i'' row in A( 5 ), the polynomial transfer matrix is proper [61]. 

6) From the structure of the input vector in (0.52), determine uft) and u^it). 
Partition polynomial matrices Eq(p) and B(p) accordingly from equation (0.63) 

Since both inputs are unknown, n^^ =0 and n^^ =2 and so E^^f s) is an empty 
matrix and Eq^{s) = Eq{s) . From the dimensions of the polynomials in equation 
(0.63), it can be seen that from the previous step, B{s) is partitioned as 



5,2(5) = 


-0.6705-0.280 2.7005 + 1.100 


0.130 


0.066 


, and 522 ( 5 ) = 


0.140 -0.640“ 
-0.054 0.240 


7) Verify By^{p) and B^^{p) are right coprime and solve Diophantine equation to 
obtain Wfp) and W^ip) in (0.64). 

In the example, 5 j 2 ( 5 ) and 522 ( 5 ) are right coprime. The solution to the 
Diophantine equation 
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E,,ip) = W,(p)B,,(p) + W,(p)B,,(p) 


are the constant polynomial matrices 


"0.000 

0.000" 


■ 0.000 

0.000" 

0.000 

0.000 


0.000 

0.000 

0.000 

-1.000 

and iTi(5) = 

-0.100 

0.064 

0.000 

0.000 


0.000 

0.000 

0.000 

0.099J 


1.500 

-0.580 


8) Form W (p) = [W^(p) WjCp)] and from (0.62) go back to the state space 
realization of the partial-state representation for simulation. 


The observer from the original system has now been successfully designed such 
that, in choosing the polynomial matrix W*{p) in (0.62), the parametrized 
Kalman Filter in polynomial form is now independent of unknown inputs MjCO • 
Changing back to state space form from the left matrix fraction realization yields 
the following system of equations: 


-0.0419 

0.9311 

0.0719 

0.0150 

-0.4610 

-0.0743 

-0.1462 

0.9586 

-0.0073 

-0.0386 

-1.0051 

-2.1930 

-1.1025 

0.0097 

0.0556 

-0.0020 

-0.2488 

-0.7881 

-0.8657 

-0.4708 

12.1667 

25.8504 

3.7543 

0.0563 

-2.3256 
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■-0.0008 

0.0000 

-0.0304 

0.0532 


0.1894 

-0.0172 

-0.0609 

0.1289 

F* = 

-0.0568 

-0.1159 

-0.3352 

-0.1925 


-0.1118 

-0.9950 

-0.0027 

-0.6358 


0.3771 

0.6721 

3.9845 

2.4990 



■-0.2405 

-0.5811 

-0.0000 

-0.9950 

0.0000 


0.1779 

2.9439 

5.1764 

-0.0449 

0.0899 


0.0000 

-0.2405 

-0.5811 

0.0000 

0.9950 


-2.5932 

-5.2463 

-0.1378 

0.0906 

-0.0016 


0.1779 

2.5882 

0.0000 

-0.0449 

-0.0000 


and ^*=0„_ 


For implementation, the resulting unknown input observer is independent of u{t) 
and the output is the full partial state (with derivatives) as 


<?(0 = ® V0 + r * j(0 

z{t) = C* q{t) + D* y{t) 


where q{t) denotes the state of the new system. 

A simulation was conducted using the above example and the estimation 
results on the original state are shown in Figure IIL7. The red line shows the 
actual state trajectory, while the other trajectories are from the three different 
observer designs. It can be seen that the polynomial UIO method (blue) outlined 
in this section performed just as good as the subspace projection UIO method 
(green). Both UIO methods, which only used output measurements to estimate 
the state, performed as well as the standard Kalman Filter, which depends on both 
inputs and outputs to estimate the state. 
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Figure III.7: State Estimate Comparison 

When deciding on a method for estimation, it is desirable to design an observer so 
that one or more inputs have no effect on the state estimation. Perhaps a system model 
has an unknown disturbance or other model uncertainties. In the case of residual 
generation, perhaps an input is assumed unknown to develop a structured residual set for 
Fault Diagnosis. When considering the subspace UIO methods and the polynomial ones 
described herein, each is designed to effectively obtain an estimate of the full state when 
the UIO rank condition is satisfied. The advantage of the polynomial methods is that 
they can still generate a residual through an estimate of the partial state. 

As mentioned before, the input replacement UIO method is limited if the 
replacement polynomials do not have higher degrees. As a consequence, this method 
works best when the degree of the numerator polynomial matrix B{p) is the same as the 
denominator polynomial matrix A{p ), which is the same as saying that D 0 in the 
equivalent state space realization. 
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The parametrization UIO method is the more general polynomial method and is 
comparable to the subspace projection UIO in complexity and performance. The main 
idea is to develop a standard Kalman Filter in polynomial form and parametrize the 
solution such that it is made independent of the unknown inputs. This method works well 
in the presence of plant and sensor disturbances with known first and second order 
statistics. One drawback common to both polynomial methods is that they require a 
linear transformation to the controllable canonical form to accommodate estimating the 
partial state. In cases where system eigenvalues are very large (i.e. greater than 10®), the 
inverse of the transformation matrix T can be near singular. 

Besides being able to generate residuals based on the partial state, another 
advantage to the parametrization method is that since the solution to the Diophantine 
equation is nonunique, it can be represented by a family of observers. If one of the inputs 
is suddenly made unavailable for the estimation, the observer can be adjusted adaptively. 
This idea is being explored as a developing possibility. 

D. INPUT OBSERVABILITY AND ESTIMATION 

Throughout the discussion so far in this chapter is the idea that if the UIO rank 
condition (0.25) is not satisfied, a partial-state estimate is still possible and can be 
obtained through polynomial UIO techniques. Specifically, it has also been shown that 
the column degrees of the denominator polynomial matrix A{p) (in the right matrix 

fraction realization H{p) = B(p)A \p)) are directly related to the ability to estimate the 
full state in (0.39). For any column degree reduction of the original A(p) to the new 
observer denominator polynomial A(p), it will result in an unknown input observer that 
cannot obtain a full-state estimate. This also implies that B(p) must have higher order 
terms. 

We have discussed residual generation based on the state and partial state. Now 
we consider it based on input and output measurements. It will be shown that in some 
cases, it is desirable to compute residuals not only on the outputs, but on the inputs as 
well. Depending on system structure, this would give added information on the class of 
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faults. It will be shown that while output estimation is straightforward when the full-state 
estimate is determined, the input estimation problem is non-trivial and it requires certain 
conditions for existence that relate to the system inversion problem. 

Compared to state estimation techniques for unknown input observers, little 
research has been conducted on input estimation [73]. Some input reconstruction 
methods base their estimation on derivatives of the output measurements [168] [169]. 
Here again, it is shown that input estimates can be determined on a larger class of 
systems from taking derivatives on the partial state. The remainder of this section 
addresses the necessary condition for input estimation and how the output equations can 
be chosen so that it can be satisfied. Similarities between the unknown input observer 
with input reconstruction and system inversion will also be addressed. Moylan [268] 
showed a method for determining the inverted system in state space. It will be shown 
here how this result relates to a system in right polynomial transfer matrix form through 
the partial state. Finally, two DC ZEDS examples will be shown to demonstrate the 
concepts. 

1. Input Observability Theory 

The goal of the input reconstruction problem is to formulate a system filter whose 
inputs are the measurements and known inputs of the original system and whose outputs 
converge to the inputs of the original system. We begin with a short example. Suppose 
we have a third order, single input, 2-output system as shown in Figure III.4 and is 
represented by the polynomial matrix fraction H(p) = B(p)A \p). If a UIO is 
developed (not necessarily polynomial-based) such that a partial-state vector estimate 
z(0 is determined, then by the polynomial matrices B(p) = [B^(p) Mp) 

defined as in (0.27) and (0.28), an estimate of the input i2(t) and output y(t) can be 
determined from (0.29)-(0.30) as 


|5^i(0= 1 ^ 0 ^ (0+ Az(t)+ Az(t)+ Az(0 

1^2(0 = 2^0^(0+ 2^1^(0+ 2^2^(0+ 
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u{t) = z{t) + a^z{t) + a^zit) + a^z{t). 


( 0 . 67 ) 


Let jiCO be arbitrarily chosen to be the “input” to the estimation u{t ), meaning that we 
solve for the highest degree term in (0.66) and substitute in (0.67) to obtain 

u{t) = ) z{t) + (^2 - ) z{t) + (<33 - ) z) + -^ y^{t), (0.68) 

1^0 1^0 


provided that is nonzero. This constraint implies that the polynomial By{p) must be 
the same order as A{p). In the state space realization, this means that the matrix D in 
the state space model must be nonzero. If is zero, then the input can be estimated 

directly from (0.67) by taking the derivative on the highest order in zXt) to obtain z (t ). 
Notice that the equation (0.68) is similar to the system inversion process, since the 
“output” is the system input estimate u(t) and the “input” is the measurement y^(t) 
[268][169]. 

Now consider the more general MIMO case, again represented in state space as 

x(t) = Ov(0 + ru(t) (0.69) 

y(t) = Cx(t) + Du(t) (0.70) 


where x{t) e , u{t) e M"", y{t) e M"" and the constant matrices O , L, C, and D are 
of appropriate dimensions. If the system is a minimal realization and a UIO is designed 
to estimate the partial-state vector, then similar to solving for the high order derivative in 
(0.66) for the SIMO case, we use the result in (0.36)-(0.37) to relate the input estimate 
u{t) in terms of the coefficients of the right matrix fraction polynomials B{p ), A{p) and 
the output y{t) as 
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am =-(bx )* ( b , - BXA)im+{ ba;' )' ym 


( 0 . 71 ) 


where e M""**"" and 6^ e are constant matrices of high degree coefficients of 

A(p) and 6(/)), respectively. Likewise, A, eM""’""' and 6^ e are the matrices of 

all lower order coefficients of A(p) and B(p), respectively. The superscript ^ 
represents the Moore-Penrose pseudoinverse. For more details on this representation, see 
Appendix D. Clearly, in order to have an inverse, A^ must be nonsingular and the 

product 5^A^' (which is the D matrix in state space) cannot be zero for a nontrivial 
solution. Again, a nonzero 5^ implies that B(p) must have the same column degree as 
A(p). If it does not, then the input estimate can be found from (0.36) as 


u(t) = A^z^(t) + A,z(t) 


(0.72) 


where z„(0 = 




P^"“ (0 is the vector of „ derivatives on the partial 


state z(t ), determined by any standard derivative filter. We have shown in Section 
IV.B.l that in cases where the UIO rank condition is not satisfied, the partial state can 
still be estimated. This result can be extended to determining an input estimate by adding 
all derivatives of all partial states up to the column degree i = I,..., n^. 


Input estimation is very closely related to system left inversion where a system’s 
inputs are determined from knowledge of its output. Necessary and sufficient conditions 
for the stable inversion problem are provided in [268]. Hou and Patton show the 
similarities between input estimation and system inversion and provide a necessary and 
sufficient condition for input observability [169]. 
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For both polynomial methods herein, special attention is paid to the input 
reconstruction problem because it provides more analytical redundancy and potential for 
improving a Fault Tolerant Control architecture. 

In order to determine conditions for input observability, consider the LTI model 
in (0.69)-(0.70) and express the system as 


-pi + O F 

x{t) 


~ 0 ~ 

C D 

u(t) 


I 


y(t)- 


(0.73) 


Now define the matrix 


M(X) = 


-n+o 

c 


F 

D 


(0.74) 


called the matrix pencil [200] [139]. Input observability is a necessary and sufficient 
condition for the existence of an estimator reconstructing inputs [169]. In the case of a 
controllable and observable system of equations (0.69)-(0.70) with unknown initial 
conditions, the input u(t) is observable if and only if 


rank(M(/l)) = dim(x)+ dim(M), VTeC. (0.75) 


If the input observability condition (0.75) is satisfied for a system, then one such input 
estimator exists of the form 


w(t) = F^w(t) + F^y(t) 


(0.76) 
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(0.77) 



where denotes the derivative of y(t) [169]. The number k is called the indev 

of the estimator and it indicates the highest derivative needed. The methods for input 
estimation described here and elsewhere require high orders of derivatives on the output 
measurements, which would make the estimate sensitive to noise. Also, in [268] an 
iterative algorithm is described for determining an input estimation filter. 

Another method for determining the input estimate is derived directly from the 
state space realization (0.69)-(0.70) and is similar to the system inversion problem in 
[268], only no iterative process is necessary. It is expressed here in partial-state space 
form as 


1(0 = r ' {(D-TEC)Tz(t) + T-'rEy(t) (0.78) 


u(t) = -ECTz(t) + Ey(t) 


(0.79) 


where all the system matrices are defined in (0.69)-(0.70), the matrix E is defined as 


E = 



if n,. = 


ifn^ >n„. 


and the nonsingular matrix T transforms state x(0 to the partial-state vector z{t) as 
before. Notice that each term in the input estimation equation (0.79) includes the matrix 
E , which is directly related to the degree of the numerator polynomial B{p) as shown in 
the equivalent equation (0.71) from the polynomial approach. Once again, in order to 
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have a proper system inverse, the degree of the numerator polynomial B{p) must be the 
same as the denominator polynomial A{p) or, in the state space realization D ^0. 

In most cases for fault diagnosis, the number of outputs will be strictly greater 
than the number of inputs. This implies that D is not square and only some of the 
eigenvalues of (^A-BEC) will depend on E. If D is rank deficient, then the inputs 

corresponding to a zero column will not have an estimate. Notice that the inverted 
“output” equation in (0.79) is identical to (0.71), which again shows the equivalence of 
minimal realizations of a system. 

If the system does not satisfy the inversion rank condition, it does not mean that 
an input estimate cannot be estimated. It only means that the system does not have a 
proper inverse. An explanation in [73] and [280] shows that if we relax the rank 
condition to 


rank(M'(/l)) < dim(x) + dim(M) (0.80) 

where 


M'(A) 


-;l/+o r 
c 0 


and we further assume that a full-state estimate is obtainable and input matrix F is full 
rank (which it will be if the system is controllable), then an input estimation process can 
be established with derivatives on the output measurements as shown in (0.76)-(0.77). 
Another way to look at this result is to consider an estimate of the input directly from the 
dynamic equation (0.32) represented in the partial state as 


u(t) = 1(0 -(t‘ r)r-‘orKo. (0.81) 
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By comparing equations (0.79) and (0.81), we see that the expression (0.81) shows 
dependence on derivatives of the partial state, but is independent of D. This relationship 
shows the similarity of determining the input estimate based on the right fraction transfer 
matrix representation in (0.72). Finally, recall that an estimate of the partial state is 
possible even when a full-state estimate is not. Then by extension, an input estimate in 
this case can still be achieved by adding derivatives as necessary to (0.81). 

2. Examples 

Two examples are presented to demonstrate input reconstruction on an unknown 
input observer. The first shows how a bank of unknown input observers successfully 
estimates the inputs while being systematically independent of each input. From the 
structured residual set, it can be seen that the indication of the presence of a fault is more 
prominent using residuals from the input estimate rather than based on the output 
measurements. In addition, for a certain class of systems the residual signal is an 
estimate of the fault signal itself. The second example demonstrates an unknown input 
observer used to estimate nonlinear components of a system. 

The first example is a UlO-based structured residual set designed to detect and 
isolate actuator faults. Suppose a system follows the flow diagram depicted in Figure 
III.8. The system model shows an input uit) being subjected to an unknown disturbance 

rj^ with known mean (zero) and variance and another completely unknown 

disturbance fit). The corrupted signal u^it) is the input to the known plant dynamics 
and its output measurement is also influenced by a zero mean disturbance with variance 
cr^. The goal is to determine the characteristics of the unknown disturbance fit) - that 

is, determine fit). 
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Figure III.8: UIO Example, Signal Flow Diagram 

To make the problem more general, suppose the system depicted in Figure III. 8 is 
a MIMO system where u{t) and y(t) are vectors. Let fit) be a scalar disturbance 
operating on one of the inputs, but it is not known a priori which one. Also let the 
disturbances rj^ and rj^, have influence on all the inputs and outputs, respectively. Now, 

the goal is to determine which input is corrupted by fit) on the basis of the input 
residual 

r^it) = uit)-uit) 

and, for comparison, the output residual 

ry{t) = yit)-yit). 

The problem is solved using unknown input observers with input reconstruction 
(UIO-IR) in a structured residual set depicted in Figure III.9 where the bank of observers 
are designed to reject the inputs one at a time. The system used in this example is the 
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sixth-order, four-input, twelve-output* port bus model for the DC ZEDS system shown in 
Figure 11.6., with parameters listed in Table II.3. For the purpose of this exercise, the bus 
is assumed to be in Mode 3 (where zones I, 2, and 3 are energized) and the system is 
adjusted in time scale using arbitrary inputs 


Mi(0 

U^{t) 

U^{t) 


sin(0 

sin(0 -I- — 
2 

sin(0 -I- n 

... n 

sm(n- 

2 . 



Figure III.9: UIO Example, Structured Residual set for Input FDI 


It turns out that twelve outputs are not necessary for the UIO implementation. A more manageable 
number is 9 (bus current, bus voltage, and capacitance voltage in each zone), which still meets all criteria 
for the minimal realization, UIO rank condtion and input observability rank condition. 
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The bank of UIOs (designed using the polynomial parametrization method) in 
Figure IIL9 shows that each of the four observers is designed to be insensitive to u^{t), 

u^{t), and u^(t), respectively, where the represents the unknown input. The 

other inputs to the observer include the outputs y (denoted as in Figure III.8) less 

the bus voltage sensors where indicated. For this problem, the system is a minimal 
realization and the UIO rank condition is satisfied for all conditions. If only the bus 
current, bus voltage, and capacitance voltage for each zone are considered as outputs, 
then the D matrix, in (0.11) would be 


D = 


0 

0 

0 

0 

0 

0 

- 0.1 

0 

0 


0 

0 

0 

0 

0 

0 

0 

- 0.1 

0 


0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

- 0.1 0 


and the input observability rank condition in (0.75) is satisfied. From this structure, it can 
be seen that ranki^D) = 3<n^ and so a solution to the input estimate can be determined 
by (0.71) or (0.79) for all except u^(t ), which shows all zeros in the fourth column. In 
other words, none of the outputs are algebraically dependent on the input u^(t ), which is 
the output voltage from the power supply in the DC ZEDS system. 

Also, since there is input dependence on the last three output measurements, they 
must be removed from the observer estimate because a fault on the input will also mean a 
fault on the outputs that are linearly dependent on the inputs. These last three rows of D 
correspond to the bus voltage sensors. 


85 



In this example, a disturbance signal f{t) was arbitrarily chosen to be a triangle 


waveform added to the input u^{t) as shown in Figure IILIO. 


Ui(t) 



Figure IILIO: UIO Example, Input Fault Trajectory 


The results from the simulation were obtained through resuiduals generated from 
the structured residual set in Figure IIL9 and show a fault indication on the far left 
column of Figures III. 11 - 111.13. The columns of these figures represent residuals 
generated from the UIOs insensitive to u^it), u^it ), u^it), and u^it) from left to right, 

where the represents the unknown input. The rows in each figure are residuals based 
on estimates of the inputs or outputs. In Figure III. 11, the subscript “ipz*” are the load 
demands leading to the SSCMs in each zone as shown in Figure II.6 as in Figure 

III.12, the subscript “ibus*” are the bus currents in each zone and in Figure III.13, the 

subscript “vb*” is bus voltage in each zone. The fourth input is from the power 
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supply but is not shown in Figure III. 11 because an estimate cannot be 
determined without derivatives to the partial state. 

It can be seen that the residual in the very top left corner of Figure III. 11 is the 
estimate of the disturbance fit) and because it resides in the first column, where no 
other signatures reside going down the column, the disturbance is on the first input uft). 
To make the figures more readable, the triangle disturbance signal fit) is superimposed 
in red, the actual results are in blue, and the mean is green. 
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Figure III.ll: Input Residuals (Fault on uft )) 


If the estimate uft) is determined from UIO.U~234 (shorthand for the unknown 
input observer insensitive to uft) ), then from the consistency between input uft ), ^ 2(0 , 
M 3 (0 , and y ^(0 , an input estimate will be uf) = uft) + f f) and the residual 
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(shorthand for input residual from UIO.U~234) would be = -f^{t), which is the 

opposite of the estimate of the disturbance signal. Any other combination of inputs and 
outputs will result in some signature in the residuals. 

Looking at the other columns in Figure III. 11, it is possible to have ambiguous 
results. For example, column U12~4 shows a similar signature in the top row with very 
small residual signatures down the column. This is essentially a problem of residual 
evaluation, but a solution might be that if the disturbance signal (i.e. the fault) is known 
to have a certain characteristics, then a matched filter can be implemented along with any 
standard noise reduction scheme to identify signals buried in noise. Another method 
might be to look for fault characteristics in the frequency domain. 
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Figure III.12: UIO Example, Output Residuals (Bus Current per Zone) 
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Looking at Figure 111.12 and Figure 111.13, we see the residuals based on output 
measurements. Once again, the very first column has no signature while all the other 
columns do (except for , which will be explained later). This means that the fault 

indication is associated with u^it). Notice in these two figures, the residuals are not 

nearly as pronounced as they were in the residuals based on input estimation. This can be 
explained by the linear combination of matrices in generating the residuals from the input 
estimation is better “tuned”, or simply more suitable, for detecting input faults. Likewise, 
residuals generated from output estimates are better suited for sensor FDl. 

As mentioned before, the only output residual with any signature in the first 
column is the one in the top left corner of Figure 111.13 (Tvi,) • This is because the bus 
voltage sensor in zone 1 is linearly dependent on the input u^it). Recall that this 

measurement was removed from U10.U~234 when determining the state estimate. This 
means that the estimate will have a u^{t) fault influence on the output measurement and 

therefore the row can be discarded. 
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Figure III.13: UIO Example, Output Residuals (Bus Voltage per Zone) 


In the second example the purpose is to demonstrate how a UIO can be used to 
estimate the nonlinear terms in a system and then further estimate unknown input. In this 
case, the power supply model with buck controller is used where its topology and 
corresponding parameters are shown in Figure II.3 and Table ILL Also, its controller 
design and parameters are shown in Figure 11.4 and Table 11.2. The power supply model 
dynamic equation is shown in (0.4) and rewritten here as 




Mi(0 
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^ 2(0 


2^ , ps 

u^{t) 


(0.82) 
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where the known inputs are \u^{t) 1 / 4 ( 0 ] = v2£' and the “unknown” inputs 

are the nonlinear terms [i/jCO u^(t)] =\d-i^ d-v^\ . The UIO design shown in 

Figure 111.14 describes the estimation procedure. The idea is to design an observer 
independent of the nonlinear dynamics, then estimate the inputs through input 
reconstruction. Finally, divide the unknown input estimates by the state estimates 

.A, ^ 

and to determine an estimate on the control input d(t). If the input can be 

measured, then a residual (t) can be formed. In this case, both the UIO rank condition 
and the input observable condition are satisfied. 



Figure III.14: Input Estimation of a Nonlinear System 


The results are shown in Figure III. 15. Here, the estimate of the control input 
d(t) is determined in the top graph while the bottom graph only serves to show that the 
estimate found is consistent with the actual control input trajectory. 
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Figure III. 15: Results from Polynomial UIO with Input Reconstruction 

From the observer point of view, if algebraic equations relating the unknown 
input to the measured output are included in the output measurements and certain rank 
conditions similar to the left inverse problem in [268], then input estimation is achievable 
without additional derivatives on the partial state. 
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IV. FAULT DIAGNOSIS APPLIED TO THE DC ZEDS SYSTEM 


It has been shown in numerous publications that the physics-based models of the 
DC ZEDS system accurately represent system dynamics [285] [335] [327] [328] [336]. 
Indeed, the DC ZEDS simulation model used herein was developed and tested by the 
authors of the aforementioned publications and the text [220]. Because the Power 
Supply, bus, SSCM, and Oring functions in Eigure II. 1 and described in Chapter II are 
well-defined, it is reasonable, then, to take a model-based approach to Eault Diagnosis on 
the DC ZEDS system. Chapter II sought to organize the DC ZEDS system into modules 
and represent their dynamics in state space form. In cases where the average value 
models are nonlinear, such as the power supply and SSCM modules, the models are 
separated into linear and nonlinear parts. Chapter III then developed several unknown 
input observer approaches that assume a model in state space form for implementation. 

The nature of the DC ZEDS system provides high bandwidth system control 
response, but because it can be controlled so tightly, faults can propagate quickly as well 
[327] (on the order of tens of milliseconds). The fault diagnosis schemes must 
distinguish between normal and casualty transients and ultimately be able to dynamically 
reconfigure under stressful battle conditions [390]. Because time is critical, model based 
methods to detect the source of a fault rather than trace its effect, is sought. 

This chapter applies a variety of model-based observer techniques (including the 
polynomial approach) to the DC ZEDS simulation model for fault detection and isolation. 
It is desired to detect faults and classify them according to their characteristics. Some of 
the means for characterization depends upon the detection method (i.e. structured residual 
set, unknown input observers, or fault modeling), while others upon the signature itself 
(i.e. positive or negative, shape, or relative size). The goal for this chapter is to develop a 
set of methodologies to identify residuals that indicate certain classes of faults. 

It is acknowledged here that standalone sensors perform a certain level of 
analytical redundancy without computations. However, individual sensors do not 
contribute to a “global” perspective as do a network of sensors connected together by 
mathematical relationships. One of the presuppositions in this work is to develop ways to 
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automate the fault detection process quickly, reliably, and accurately - and a hierarchical 
automated fault detection layer using robust model-based estimation is one such way. 


A. MODULAR APPROACH TO INPUT-OUTPUT CONSISTENCY 

The DC ZEDS model as a whole contains 68 states, and this does not consider the 
ac side of the overall power system from prime mover to constant power load. From the 
input of the power supply to the input to the SSIM in each zone, the DC ZEDS system 
has many modes and system configurations that make it difficult to estimate as a whole - 
even under normal operating conditions. A suitable solution to handle the sheer size of 
the problem is to separate the system into smaller logical subsystems and connect them 
together with common inputs and outputs. This modular approach to verifying input- 
output consistency is done for the DC ZEDS system and the topology is shown in Figure 
IV. 1. 



Figure IV.I: DC ZEDS Module Topology 
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Four representative component module block types are configured so that each 
input/output pair between them is shared measurements such that the output to one block 
is the input to the connecting block. Under this method, the entire system is viewed as a 
dynamic interconnection of subsystems with common inputs and outputs and is 
consistent with the DC ZEDS simulation model. Notice also in Figure IV. 1 that the 
models show voltages as outputs moving toward the constant power load in each zone 
and currents as inputs. 

Using various types of observers, model consistency is verified by comparing 
module estimates with measurements. With algebraic equations relating the states and 
the input in the form of energy conservation laws, the observers become input observable 
and therefore an estimate of the inputs can be made without derivatives. Taking a closer 
look at the input/output connections shown in Figure IV.2, we see that redundancy is 
verified from the input and output measurements of each module block “A”, “B”, and 
“C”. 



Figure IV.2: System Module Interconnection 


Suppose block “B” is the port SSCM in Zone 1 where the external inputs and 
in (0.19)-(0.20) are the port bus voltage and load current and let and be the 

SSCM output voltage and the Zone 1 load requirement to the port bus, respectively. To 

verify analytical redundancy across the entire DC ZEDS system, estimates can be made 
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on system measurements from both sides of the block and compared to the estimates (or 
measurements) from each connecting block. In our example, and from the port 

bus (block “A”) would represent the bus voltage and load current in (O.lO)-(O.ll), while 
and Je, from the Oring function (block “C”) would represent input and output 

^pon (0.21)-(0.22). Residuals can be generated by comparing the measurements from 
each common line as 


rBc(t) = \(t)-ybS.t) 
resit) = It hSt)-ycSt) 


Notice that the one estimate in each equation of (0.83) is independent of the other. 
In this manner, the modular interconnection model shown in Figure IV.2 can zoom down 
to any desired level in the system, provided its external inputs and outputs are the outputs 
and inputs to connecting modules. 

Among the benefits of this consistency scheme is increased redundancy. For 
example, if is a control signal input for the controller in Module A, and the signal 

becomes disconnected or otherwise erroneous, then between the two estimates u„ and 

“2 

the fault can be isolated to one of the modules. In some cases, different methods of 

FDI can be performed in separate modules that draw the same conclusion and your 
confidence of the fault indication becomes much stronger. 

From Module A, verifying sensor consistency is done using a structured residual 
set of regular Kalman Filters as shown in Figure IV.3. Here, an observer is designed to 
minimize the mean squared error of the output residuals using all the inputs and all the 
outputs less one. Each successive observer in the set performs the same operations until 
all Hy sensors are removed one at a time for a total of observers. The notation ^y(t) 
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indicates the output measurement vector y(t) without sensor 1, and so on. Likewise, 
^y(t) and are the corresponding output measurements and output residuals based 

on observer inputs u(t) and y(t) without sensor 1. The residual generator compares the 
full output vector y(t) with all its estimates for consistency. Suppose for example that 
the number three output sensor is faulted. 


.v(d 



Figure IV.3: Sensor FDI Structured Residual Set 

Then by the observer structure in Figure IV.3, the third observer in the set would 
accurately estimate the faulted output signal 3, while the rest would show a signature in 
the residuals. This is because all the other observers are dependent on the faulted 
measurement. 
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The following DC ZEDS example demonstrates the concept of modular EDI. 
Consider the nominal operating scenario where the port and starboard power supplies are 
energizing the port and starboard bus, respectively. The starboard bus is providing power 
to the load in Zone 1 while the port bus is feeding the load in Zones 2 and 3 with 
operating parameters according to Table 11.5. The power supply controller is regulating a 
500 V output voltage through the control input d , which is assumed to be measurable. 

Suppose a fault occurs to the load current signal in the port power supply 
controller (shown in Eigure II.4) at time = 0.6 seconds. In the fault condition the 
signal is a constant zero, thus removing the feedforward path to the PI controller. 

The local effect of the fault is that the feedforward path will communicate to the 
controller that there is no required load current, which disables the controller’s ability to 
provide fast response to changes in the load. The fault propagates through to the port bus 
as a small transient that quickly recovers to the steady state as shown in Eigure IV.4 and 
Eigure IV.5. The starboard bus remains unchanged. 
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Figure IV.4: Bus Voltage Showing Bus Current Sensor Fault 
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Starboard Bus: Bus Current, Zone 1 



Figure IV.5: Bus Current Showing Bus Current Sensor Fault 


The goal here is to use the modular FDI approach to detect the fault and 
determine in what module the fault originated. Beyond that, it is desirable to gain other 
fault characteristics, if possible. Looking at the power supply module shown in Figure 
IV. 1 we see that it is connected to the port bus through common signals and . 

The first step to our goal is to conduct sensor FDI on the sensors within the bus module. 
Developing a structured residual set of Kalman Filters as shown in Figure IV.3 achieves 
that task. So as to not be unnecessarily exhaustive in showing all the sensor FDI 
residuals, only the zone current sensors are plotted here since the other residuals are 
consistent with what is shown here in Figure IV.6. 
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Figure IV.6: Bus Sensor FDI Residuals for Bus Current Sensor Fault in Power 

Supply 

The residuals in Figure IV.6 show that there is consistency between estimates and 
measurements since they are white and stationary. It can be shown that the other sensors 
in the set show the same result. Based on this result, we conclude that the sensor fault is 
not in the bus module. Now we take a look at the structured residual set in the port power 
supply module. Since the fault in question occurs on the input to the power supply, an 
actuator fault diagnosis problem similar to the one described in Figure III.8 applies. 
Applying the structured residual set similar to the one shown in Figure III. 9 for the power 
supply model in (0.4)-(0.5), we receive the following residuals on the input and output 
estimates as shown in Figures IV.7 and IV.8. 
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Figure IV.7: Input Residuals for Bus Current Sensor Fault in Power Supply 

Controller 


In contrast to the sensor residuals in the bus module in Figure IV.6, the residuals 
in Figure IV.7 indicate the presence of a fault. As before, the residuals are generated 
from the input estimations indicated by row and the unknown input observers by column. 
Looking down column U123~, it can be seen that the only residual which is not white is 
in the last row, which corresponds to the signal in question, . This makes sense 
because the estimates in the column to the far right base their estimate on all inputs 
except , which is where the fault resides. The residuals down the column indicate 

input/output consistency, while the input estimate for contains the fault signature. 
Similar results are shown in the output residuals in Figure IV.8. 
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Figure IV.8: Output Residuals for Bus Current Sensor Fault in Power Supply 

Controller 


The result of the bus current sensor fault example shows that first, the fault did 
not have any indication of originating in the bus module. It did, however, have an 
indication in the power supply model and, more specifically, the fault showed indication 
of residing in the bus current sensor somewhere between the bus and the controller. 

B. FBI BASED ON FAULT MODELING 

In the last section we considered a class of faults where the fault information was 
due to a difference between measurements and estimates. Using the methodical, 
structured residual set for both input and output residuals, a fault was detected. The key 
component to this method is a valid model that accurately represents the input-output 
relationship. In this section, we consider another means to detect faults using 
mathematical models specifically “tuned” to a particular class of fault. In this case. 
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regular Kalman Filters are determined from models which are known to indicate a fault 
condition in the DC ZEDS system. In this case, indication of a fault in the residuals is 
exactly opposite of the results from the previous section. A zero residual would indicate 
a match to the fault model, providing precise fault detection and isolation - provided the 
fault model is consistent with the actual fault dynamics. 

Suppose we have a nominal plant configuration similar to the one described in the 
last section and we wish to consider what the effects of an SSCM switch failure in the 
closed position would do to the DC ZEDS system. Then, the goal is to model this failure 
and simulate the EDI scheme. 

Eor simulating the SSCM switch in the closed position, we first choose the 
starboard SSCM in Zone 1 depicted in Eigure ILL In the average value model, the fault 
is triggered by setting the SSCM control signal d to 1. Similarly, in the detailed model, 
the control switching signal s is set to 1. A simulation of the fault condition was 
conducted and its effect on the bus is shown in Eigures IV.9 and IV. 10. 
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Figure IV.9: SSCM Switch Closed Failure: Port/Stbd Bus Voltage 


From Figure IV.9 it can be seen that the SSCM failure has very little effect on the 
bus voltage. This is due to the sophisticated and robust power converter control elements 
that provide nearly-ideal load regulation and transient performance [327] [332]. The port 
and starboard bus current is shown in Figure IV. 10 where the starboard bus shows a large 
transient at the time of the fault, but settles back to near the same level as before. It can 
also be seen on the starboard bus current figure that the average value model follows very 
closely to the detailed model, even in the fault transient. Notice also that the port bus 
remains unchanged. 
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Figure IV.IO: SSCM Switch Closed Failure: Port/Stbd Bus Current 


The next set of figures shows the inputs and outputs to the SSCM under the 
switch closed fault condition. Figure IV. 11 shows the inputs to the SSCM. The bus 
voltage is virtually unchanged as before, but the load current demand oscillates very 
rapidly in the steady state. This is due to the constant power nature of the load that the 
SSCM is feeding. If the voltage to the load oscillates in the steady state, then so will the 
current according to the relationship = P *, where P* is the constant load 

requirement shown in Figure 11.15. Here again, the average value fault simulation agrees 
with the detailed simulation. 
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Figure IV.ll: SSCM Switch Closed Failure: SSCM Inputs 


Looking at the SSCM circuit diagram in Figure 11.10, it can be seen that if the 
switch is closed permanently, then the steady-state output should be very close to the 
input. This is indeed the case as shown in Figure IV. 12 where the dc value of the SSCM 
output is close to the bus voltage (500 V). 
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Figure IV.12: SSCM Switch Closed Failure: SSCM Outputs 

Next we examine what is happening inside the controller depicted in Figure 11.12. 
As it turns out, the position of the control signal d becomes important for detecting the 
fault. First, it makes sense that the duty cycle d is around 80% under no-fault conditions 
because the regulated output voltage is 80% of the input voltage. 

If the SSCM switch fails closed, then we see that the output voltage becomes 
close to the input voltage. The controller in Figure 11.12 shows a large error signal 
between the output voltage (-500 V) and the reference voltage (400 V), resulting in the 
controller commanding the switch to open. This fault condition means that the command 
signal J is a constant zero, while the actual operating condition of the switch is stuck in 
the opposite state. These command signal dynamics are shown in Figure IV. 13. 
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Figure IV.13: SSCM Switch Closed Failure: SSCM Control Signal d 


In what follows we derive the equations that model the fault. Consider the SSCM 
model equations in (0.19)-(0.20) where the dynamic equation is rewritten here as 


i = A, 


5, 


M, 


Ur. 

1 

2,sscm 

L 

_U^_ 



(0.84) 


Recall that the inputs u^it) and u.^{t) are the nonlinear terms relating the control signal 
d and the states of the system as 


u^it) = d ■x^{t) 
Uj(t) = d ■x^{t) 

109 


(0.85) 




































Under the SSCM switch closed fault condition, the control signal d = l and the 
nonlinear terms in (0.84) can now be expressed as a part of the linear dynamics in the 
new fault dynamics state space equation 


^( 0 = 


l,sscm 


Ml(0 

U^(t) 


( 0 . 86 ) 


where 


0 0 0 0 
0 0 0 0 


Al^ = A„, 


+ 


0 0 0 0 


0 0 0 0 

0 0 — 0 

u... 


0 

0 

1 

0 

0 


The fault model in (0.86) along with the output equations in (0.20) is controllable, 
observable, and input observable. Therefore, a standard Kalman Filter is applied to the 
SSCM switch closed fault model and simulated using the detailed model data. The 
results in Figure IV. 14 show the opposite trend in a positive identification of the fault 
compared to previous methods, which is expected. What we have is essentially a filter 
that is “tuned” to the particular fault where the residuals will be zero when the fault 
occurs, which indicates input/output consistency in the fault condition. All five states of 
the SSCM model reveal a constant near-zero mean residual, but only 4 of the 5 are shown 
in Figure IV. 14. 
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Figure IV.14: SSCM Switch Closed Failure: Fault Model Residuals 


Finally, we go back to the nonlinear input estimation problem described in Figure 
111.14. The SSCM model (0.19) has a similar model equation structure as the power 
supply model (0.4) and so its control signal d can be estimated in like manner. Recall 
that in the failed condition where the output voltage is far greater than the reference 
voltage, the SSCM controller is continuously sending the command to have the switch 
opened as seen in Figure IV. 13. Without the use of the fault model, a fault indication can 
still be achieved if a UIO scheme is set up such that it estimates the command signal d . 
If the control signal is measurable, then a residual on the command signal can indicate the 
presence of an SSCM switching fault as shown in Figure IV. 15. 
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SSCM Switch Fault (Closed): Residual on SSCM Control Signal D 



Figure IV.15: SSCM Switch Closed Failure: Residual on Control Signal d 

One more example for FDI based on fault modeling is presented. Suppose only 
the port power supply is on line providing power to all three zones and an open circuit on 
the bus occurs so that Zone 3 is dead. The scenario was simulated and the results are 
shown in Figure IV. 16. The bus current, is shown to drop from the no fault Mode 3 
condition down to the Mode 2 condition. This is consistent with what is expected since 
the load in Zone 3 is completely dropped. 
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Figure IV.16: Bus Open Circuit in Zone 3: Bus Current 


We use the three bus mode models shown in Figure II.6 and Table 11.4 as a means 
to detect the bus open circuit. Using the inputs and output measurements, we process the 
data through three Kalman Filters tuned to each bus mode configuration. The model that 
is most consistent with the input and output data will show a zero residual. The observers 
are based on bus model equations (O.lO)-(O.ll), and the results are shown in Figure 
IV. 17. From the residuals, it can be seen that the system began in Mode 3 operation (i.e., 
all zones on line) and then suddenly switched to zone 2, indicating an open circuit in zone 
3. 
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Figure IV.17: Bus Open Circuit: 3 Mode Residuals 


To conclude this section, it is shown that FDI using fault models shows good 
results. But like all model-based methods, they are dependent on the accuracy of the 
model. Fortunately, these models have been proven to be accurate. 
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V. CONCLUSION 


A. SUMMARY OF FAULT DIAGNOSIS RESULTS TO THE DC ZEDS 

MODEL 

Model-based Fault Diagnosis is an excellent method for validating system 
input/output consistency and to detect fault indicators in the residuals, which is the 
fundamental principle of analytical redundancy. However, these methods are only as 
good as the models upon which they are based. Fortunately, it has been shown that the 
DC ZEDS model is an accurate power electronics-based distribution system and is very 
similar in design to the reduced scale hardware model implemented in the Naval Combat 
Survivability (NCS) Testbed located at Purdue University, which is specifically designed 
for concept studies such as the ones in this work [335]. The DC ZEDS system is a 
reduced-scale model, yet it contains enough of the complex, nonlinear, dynamically 
interdependent nature to make it quite suitable for scholarly research 
[56] [327] [328] [245] [246] [335] [285]. 

Two models were used side-by-side for the results of this research: the average 
value model and the more accurate detailed model. The difference between the two 
models is the 20 kHz switching dynamics of the power supply and SSCM. The average 
value models are nonlinear, but the switching dynamics are removed and, through the use 
of unknown input observers, accurate state and input estimates are determined for robust 
residual generation. In these fault scenarios, most of the EDI results aim at detecting the 
source of a fault rather than studying its propagation effects. 

Based on the results, the average value model accurately matched the detailed 
model for the fault scenarios studied herein - even in the transients. It is acknowledged, 
however, that a class of faults caused by switching dynamics may not be accurately 
modeled (and therefore not accurately estimated) using the average value model and so 
they were not addressed. 

The closest one of the fault scenarios came to studying switching dynamic faults 
is the one with the SSCM switch failure in the closed position. Here, the nonlinear 
dynamics are essentially removed by a constant connection, and a model of the fault is 
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formed. The FDI scheme in this scenario used an observer based on the fault model to 
generate a residual. Similarly, another fault model scheme presented in the results uses a 
bank of Kalman Filters specifically “tuned” to detect an open circuit in various places 
along the bus. Again, the indication of this fault will be a zero in the residuals, meaning 
that there is input/output consistency in the fault condition. As long as the models 
accurately represent the system and the fault, the fault model FDI scheme works quite 
well. One drawback to this method might be that it is specifically designed to detect a 
particular fault while other FDI methods (using a structured residual set, for example) are 
able to detect a broader class of faults. This might not be undesirable, though, if several 
FDI methods are applied concurrently; for example, one FDI method might be used for 
detecting a broader class of faults with less accuracy and the other method might be used 
for detecting (with greater accuracy) a very specific fault that is known to produce grave 
consequences. 

For a more general FDI scheme, it is shown that model consistency for the overall 
large-scale system can be successfully subdivided into modules with common inputs to 
and outputs from connecting modules. In this manner, classes of faults (such as sensor 
faults or communication faults between modules) can be isolated to a particular module 
by verifying input/output consistency between modules. This modular FDI approach has 
the effect of being able to subdivide into arbitrarily small submodules, provided that the 
inputs and outputs are common to the connecting modules. 

B. SUMMARY OF POLYNOMIAL UIO RESULTS 

A special type of observer that can be made insensitive to any combination of 

inputs (provided certain conditions are satisfied) is used extensively for the FDI 

application. It is shown in several examples that a structured residual set of unknown 

input observers can accurately detect and isolate actuator faults and in some cases can 

estimate the fault signal itself. Here, two UIO methods based on a polynomial 

representation of a system are developed and compared to the subspace projection UIO 

method. Since these methods take advantage of the partial-state representation of a 

system, it is shown that in some cases an estimate of the partial state can still be 

determined even if the UIO rank condition is not satisfied, which is a necessary condition 
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for the subspace projection UIO method. For the purpose of residual generation, a 
partial-state representation could possibly be all that is needed. If the full-state estimate 
is required, then it can still be done by taking derivatives of the partial state - a process 
that is sensitive to noise. 

The input replacement UIO method is straightforward to implement, but is 
restricted to a smaller class of systems where higher degrees of the numerator polynomial 
is necessary. One way to ensure the higher numerator polynomial degree is to include 
output measurements that are dependent on the states and inputs. In the state space 
realization, this implies that the D matrix is nonzero. The input replacement UI method 
is also similar to the input inversion problem. 

The parametrization UIO method is comparable in performance to the subspace 
projection method and applies to a larger class of systems. It is based on the nonunique 
solution of the Diophantine equation, where the particular solution we choose separates 
the known system dynamics from the unknown. It is believed that this concept can be 
extended to an adaptive polynomial UIO method where the parameter vector is 
determined adaptively by minimizing the output residual. However, precise conditions 
for existence of a solution still have to be determined. In some applications to the DC 
ZEDS models, the parametrized solution of the UIO provided a nonrealizable estimator. 

Having a feed-through D matrix in the system of equations of the state space 
realization (which implies higher order numerator polynomials in the right matrix 
polynomial realization) also allows for input estimation without having to take 
derivatives of the partial state. Input estimation is quite useful for fault diagnosis as 
shown in several examples where the input to a nonlinear model was estimated. It is also 
a key element to the modular approach to EDI as it enables each module to estimate both 
outputs and inputs for consistency validation. It has also been shown that in some cases 
an estimate of the fault signal itself can be determined. 
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C. FUTURE WORK 

For follow-on work, other non-model based fault diagnosis strategies could be 
developed and compared to these results. It would also be valuable to validate the fault 
diagnosis schemes using the Naval Combat Survivability (NCS) Testbed; it is very 
similar in design to the DC ZEDS model described in Chapter II and has been developed 
by the same team of engineers at Purdue University [335]. Here, real-time EDI on an 
actual hardware DC ZEDS system could be compared with simulation results. The work 
would require modeling and simulation on a validated truth model using Advanced 
Continuous Simulation Eanguage (ACSE). 

At the reconfiguration level shown in Eigure I.l, raw residual data is analyzed and 
decisions are made about positive fault identification. Erom the innovations determined 
in Chapter IV on the stochastic models described in Chapter II, residual evaluation 
techniques can be applied to determine fault probabilities or to develop an intelligent 
fault decision-making process. Methods for residual evaluation may include: 


• Generalized likelihood ratio (GER) testing 

• Euzzy logic 

• Adaptive thresholds (time and/or freq domain) 

• Hypothesis testing 

• Multiple Model Adaptive Estimation (MMAE) 


In this systematic manner of evaluating residuals, a positive fault indication can then be 
sent to a control reconfiguration scheme to achieve the ultimate goal: continuity of 
service despite combat damage involving cascading failures. 

Other areas for EDI research could be to introduce multiple simultaneous faults or 
conduct similar fault simulations on the DC ZEDS system with a droop a mode capability 
where power in each zone is brokered between the port and starboard busses. 
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In [60], authors Chen and Patton note that robust FDI is still an “open problem” 
where the most significant challenge is due to conflicting criteria in model-based robust 
residual generation. Specifically, residual generation methods must be sensitive to faults, 
yet insensitive to modeling uncertainties. In order to account for this in the residuals, a 
multicriteria optimization approach can be investigated where fault effects are maximized 
and modeling uncertainties are minimized. Some general optimal FDI work has been 
studied using a single objective function codified as a weighted sum of objectives 
[59] [60]. In many of these cases, however, the solutions do not adequately address the 
original problem [322]. What’s more, FDI is most often considered separately from the 
rest of the system design. A better approach might be to develop a multicriteria optimal 
solution by investigating trade offs between design variables, system constraints, control 
laws, sensor allocation, and fault sensitivity. 
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APPENDIX A: IFAC SAFEPROCESS FAULT DIAGNOSIS 

DEFINITIONS 


The following keywords are common terms applied to the field of Fault Diagnosis and 
Fault Tolerant Control. They were developed by the Fault Detection, Supervision, and 
Safety for Technical Processes (SAFEPROCESS) Technical Committee and published in 
Isermann and Bade, 1997; [178]. The definitions are widely accepted in published books 
on the subject (Blanke et. al., 2003; [44], Chen and Patton, 1997; [60]). 


Active fault-tolerant A fault-tolerant system where faults are explicitly detected and 
control system accommodated. Opposite of a passive fault-tolerant system. 


Analytical 

redundancy 


Use of two or more, but not necessarily identical ways to determine a 
variable where one way uses a mathematical process model in analytical 
form. 


Availability 


Constraint 


Dependability 


Diagnostic model 


Probability of that a system or equipment will operate satisfactorily and 
effectively at any point in time. 

The limitation imposed by nature (physical laws) or man. It permits the 
variables to take certain values in the variable space. 

A form of availability that has the property of always being available when 
required. It is the degree to which a system is operable and capable of 
performing its required function at any randomly chosen time during its 
specified operating time, provided that the item is available at the start of 
that period. 

A set of static or dynamic relations which link specific input variables - the 
symptoms - to specific output variables - the faults. 


Discrepancy 


An abnormal behavior of a physical value or inconsistency between more 
physical values and the relationship between them. 


Disturbance 


An unknown (and uncontrolled) input acting on a system 


Error Deviation between a measured or computed value (of an output variable) 

and the true, specified, or theoretically correct value. 

Fail safe The ability to sustain a failure and retain the capability to make a safe close¬ 

down. A system where the occurrence of a single fault can be determined 
but not isolated and where the fault cannot be accommodated to continue 
operation. 
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Fail-operational 

Failure 

Failure effect 

Failure mode 
Fault 

Fault 

accommodation 
Fault detection 
Fault detector 
Fault diagnosis 

Fault estimation 
Fault 

identification 
Fault isolation 

Fault modeling 

Fault propagation 
analysis 

Fault tolerant 
control 


A fault-tolerant system has the property that a single fault does not develop 
into a failure of the closed-loop system. A Fail-operational system is one 
that the performance of the faulted system remains the same (original 
objective of the nominal system is met). 

Permanent interruption of a systems ability to perform a required function 
under specified operating conditions. 

The consequence of a failure mode on the operation, function, or status of an 
item. 

Particular way in which a failure can occur. 

Unpermitted deviation of at least one characteristic property or parameter of 
a system from its acceptable/usual/standard condition.[Blanke et. al. adds: 

A fault is the occurrence of a failure mode (Blanke, 2003)] 

A correcting action (reconfiguration or a change in the operation of a 
system) that prevents a certain fault to propagate into an undesired end- 
effect. 

Determination of faults present in a system and time of detection. 

An algorithm that performs fault detection and isolation. 

Determination of kind, size, location, and time of occurrence of a fault. 

Fault diagnosis includes fault detection, isolation, and estimation. 


Determination of the size and behavior of a fault over time. 

Determination of the size and time-variant behavior of a fault. Follows fault 
isolation 

Determination of kind, location, and time of detection of a fault. Follows 
fault detection. 

Determination of a mathematical model to describe specific fault effect. 


Analysis to determine how certain fault effects propagate through the 
considered system. 

Aims at changing the control law so as to cancel the effects of the faults or 
to attenuate them to an acceptable level. 


122 



Fault tolerant 
system 

Hardware 

redundancy 

Incipient fault 

Malfunction 

Monitoring 


Passive fault- 
tolerant system 

Perturbation 


Protection 


Qualitative 

model 

Quantitative 

model 

Recoverability 

Reliability 

Remedial action 


A system where a fault is accommodated with or without performance 
degradation, but a single fault does not develop into a failure on subsystem 
or system level. Fault tolerant system has the property that a single fault 
does not develop into a failure of the closed-loop system (p. 9). 

Use of more than one independent instrument to accomplish a given 
function 

A fault where the effect develops slowly e.g. clogging of a valve. In 
opposite to an abrupt fault. 

An intermittent irregularity in the fulfillment of a system's desired function. 


A continuous real-time task of determining the conditions of a physical 
system, by recording information, recognizing and indicating anomalies in 
the behavior 

A fault-tolerant system where faults are not explicitly detected and 
accommodated, but the controller is designed to be insensitive to a certain 
restricted set of faults. Contrary to an active fault-tolerant system. 

An input acting on a system, which results in a temporary departure from the 
current state. 

Means by which a potentially dangerous behavior of the system is 
suppressed if possible or means by which the consequences of a dangerous 
behavior are avoided. 

Use of static and dynamic relations among system variables and parameters 
in order to describe a system's behavior in qualitative terms such as 
causalities or if-then rules. 

Use of static and dynamic relations among system variables and parameters 
in order to describe a system's behavior in quantitative mathematical terms. 

Possibility to accommodate the fault or to reconfigure the system if fault 
occurs. 

Ability of a system to perform a required function under stated conditions, 
within a given scope, during a given period of time. 

A correcting action (reconfiguration or a change in the operation of a 
system) that prevents a certain fault to propagate into an undesired end- 
effect. Synonymous to fault accommodation. 
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Residual 

Safety 

Safety system 

Sensor fusion 
Supervision 

Symptom 

Threshold 


A fault indicator, based on a deviation between measurements and model- 
equation-based computations 

Ability of a system not to cause danger to persons or equipment or the 
environment 

Electronic system that protects local subsystems from permanent damage or 
damage to environment when potential dangerous events occur. 


Integration of correlated signals from different [sources]. 

Monitoring of a physical system and taking appropriate actions to maintain 
the operation in the case of faults. 

A change of an observable quantity from normal behavior. 

Limit value of a residual's deviation from zero, so if exceeded, a fault is 
declared as detected. 
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APPENDIX B: POWER SUPPLY AVM TRANSEORMATION 


The average value model begins with several reference frame transformations of 
the input voltage to obtain the time-dependent rms voltage of the input. The source for 
the transformation is shown in Krause et al [220]. 

The 3-phase ac input voltage is expressed as 


ag 


= V2£'cos(6'^^ 


= 'JlE cos 




eg 


= 42 E 


f 


cos 


vf 


(.87) 

(. 88 ) 

(.89) 


and the arbitrarily rotating real and reactive (q-d) axis representation is expressed as 


^gdOs ^abc ’ 


(.90) 


where 


cos 6*^ 

cos(^^ 

-Inji) 

cos (6*^ +2nl3) 

sin^^ 

sin(^^- 

-In 13) 

sin (6*^ +2nl3) 

1/2 

1/2 

1/2 


If the input deviates from being balanced, then the third of 3 outputs of the 
transformation will be nonzero. 
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A second transformation is applied to the arbitrarily rotating orthogonal q-d axis 
to a reference frame where the reactive voltage is identically equal to zero. In doing so, 
the transformation calculates the rms voltage of the input, yflE. 




V2E 



0 


I ^g^qdOs 


(.91) 


where 


a 




cos(^^) Sin(^3^) 0 

-sin(^3^) cos(^3^) 0 

0 0 1 


and 


g’g = tan 




(.92) 


The average value model of the 3-phase bridge rectifier is developed in [Krause, 
et al] and the schematic of the Port Power Supply in the DC ZEDS model is shown in 
Figure 1 along with its Buck Converter controller. The input voltage expression takes 
into account the line-to-line voltage of the balanced 3-phase ac input and represents it 
with the rms voltage. The 3 phase full wave rectifier equipped with thyristors allows for 
full control of the output dc voltage through the control input firing angle a . 
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APPENDIX C: SUBSPACE PROJECTION UNKNOWN INPUT 

OBSERVER 


The purpose of Appendix C is to develop the subspace method Hui and Zak’s 
unknown input observer from a slightly different approach presented in [175]. The 
resulting observer design allows for any combination of inputs to be made insensitive to 
the observer - or all of them - provided that certain rank conditions are met. This design 
freedom makes this observer type useful for Fault Diagnosis. The paper does not 
explicitly address output measurements that are dependent on the inputs (i.e. D^O), 
only to the extent that the input can be subtracted from the output measurements before 
the UIO is implemented and then added back to obtain an output estimate. 

The class of dynamical models considered herein is of the form 


X = Ax + Bu 

(.93) 

II 

(.94) 


where state vector x e M"', input u e M"", and output y e M"" and matrices A , B and C 
are completely known with appropriate dimensions. For a system where at least one of 
inputs is unknown, the input matrix B and corresponding input vector function u are 
partitioned into known and unknown inputs as 


Bu = [5j 



where the first known inputs are separated from the remaining unknown 

inputs. In this manner, the inputs of ( .93) can be represented with separated known and 
unknown inputs with the same output equations as 
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(.95) 


x = Ax + B^Uy + 
y = Cx 


where B^ e M""‘ and B^ e M"“" . It should be noted up front that the well-known rank 
condition must be satisfied in order to check for the existence of the observer [230]. For 
the observer to exist, rank(C52) = rank(52), which implies that there must be at least as 
many independent outputs as there are unknown inputs. 

The Hui and Zak UIO algorithm is based on the idea that we want to project the 
unknown component in M"' onto the subspace of known state components. Then, an 
observer is constructed within this subspace to estimate the entire state. Let F be a 
projection matrix such that 


P = MC 


(.96) 


where P e and M e . The projection matrix has two important properties. First, 
P is symmetric {P^ =P)- Second, the square of P is itself (P^ =P)- The second 
property also implies that 


P" = P Vn > 1 


since any positive integer exponent n can either be divided by two with a remainder of 0 
or 1. When expanding the exponents of P the result will always be a series of 
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multiplications of the projection matrix raised to the first or second power, which will 
always reduce to the second property of the projection matrix. Using the projection 
matrix, the state vector can be separated as 


x = {^l - P^x + Px 


(.97) 


Notice that the expression (^I -P^ is a projection matrix itself because it is symmetric 
and 


[I-Pf -IP + P^ 
=I-2P+P 


=(/-/>). 


We can also conclude that 


P{I-P) = P-P^ =0. 


Substituting the output equation (.94) and (.96) into equation (.97) we get 


x = q +My 


(.98) 


where 


q = {l-MC)x. 


(.99) 
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The projection matrix I -MC is constant so the derivative of (.99) is simply 


q = {l-MC)x 


(.100) 


Notice that q can be viewed as the projection (not necessarily orthogonal) of x onto the 
subspace of known components in M"' . If the subspace state q can be estimated, then an 
estimate of the entire state x can be determined by equation ( .98) since the output 
measurement vector y is known. To determine M , the constraint 


{1-MC)B^=0 


(. 101 ) 


must be met, which essentially removes the influence of the unknown input in ( .95). An 
observer is then designed within the subspace. Provided that the rank condition 
rank{CB^^ = ranki^B^^ is met (Kudva et al [230]), the known state subspace exists and 
the matrix M is determined to be 


M=52((C5,)'+ //„(/,-(C5,)(C52)')) (.102) 


where the superscript t denotes the Moore-Penrose pseudo-inverse operation and 
Hq is ^ design parameter matrix introduced in [175] (default is zero). To 

determine the dynamic state space model for the subspace state q , we substitute equation 
(.95) into (.100) to get 

q = {l-MC)A[q + My) + [l-MC)B^u,+{l-MC)B^u^ ( .103) 
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Recall that M is chosen such that the constraint (.101) is met and so the last term 
in equation - the term containing all unknown inputs - is zero. The resulting state space 
equation is of the form 


q = Aq + Bu (.104) 

y = Cq (.105) 


where 


A = {I-MC)A 
B = {I-MC)[AM fij] 


u = 


y 

u. 


and 


j j - CMy 


Notice the new input u is the vector function of output measurements and known 
inputs, which is fully known. From here, it can be shown that a standard proportional 

gain observer can be applied to equations (.104)-( .105) provided the pair is fully 

observable. Note that it is possible that the original system might be fully observable 
while the pair C j is not. Choosing the closed loop gain matrix L to be 


L = {I-MC)L, 


where L is the observer gain matrix chosen for the original system in ( .93) and ( .94). 
Finally, the estimate of the subspace state q is determined to be 
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Now to show convergence of the subspace method, choose the state error to be e = q-q 
so that the error dynamics are determined by the equation 


e = {A — LC\e. 


Knowing that A = (^I -MC)A and L = (^I -MC)L defined in equations ( .104) and 

( .105), the error dynamics can be expressed as the error of the original observer 
projected onto the known system subspace as 


e = {l-MC){A-LC)e. 


(.107) 


Looking at equation ( .107), it is not readily apparent that the error converges to zero. 
With a few observations, though, it can be shown that it does. First, pre-multiply 
equation (.107) by the projection matrix MC , which yields 


MCeit) = 0 Vt 


(.108) 


since 


MC{l-MC) = MC-{MCf =MC-MC = 0, 



which implies that 


MCe(t) = constant Vt. 


Next, pre-multiply the projection matrix MC on q in equation (.99), yielding 


MCq{t) = 0 Vt 


for the same reason as in ( .108). Choose the initial condition on the estimate q in the 
expression as 


MCq{0) = 0, 


then 


MCe{t) = 0 Vt. 


For simplicity, define a projection matrix P' and system matrix A! as 


P' = I-MC and 
A' = A-LC 


then the expression for the error equation in (.107) can be expanded as 
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f 


P'e(t) = P' 


I + P'A't + 


- + ... 


2 ! 


>( 0 ) 


= P' 


r V ^ 

I + At + -+ ... 


V 


2 ! 


.( 0 ) 


J 


= P'e^'‘e(0) 


Now choose L so that the desired plant characteristics A' = (A-LC) converge to 
zero so that (/ -MC)e(t) 0 . Putting it all together, it can be seen that the error e(t) 
approaches zero because 


{I-MC)e(t) = e(t)-MCeit) 0 
^ e{t) 0. 

Finally, since it has been shown that the estimate q of the known component 
subspace state converges to the actual state q , an estimate x for the whole state can be 
determined based on the relationship between q and v in equation (.98) as 


x = q + My. 
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APPENDIX D: POLYNOMIAL MATRIX DELINITIONS AND 

THEOREMS 


The purpose of Appendix D is to provide definitions to polynomial matrices for 
which the multivariable unknown input observer method is based. Most of the 
definitions herein are documented in [61],[200], and [229]. 


Definitions 

Definition D.l (field). Afield consists of a set, denoted by T, of elements called scalars 
and two operations called addition " +" and multiplication " •"; the two operations are 
defined over T such that they satisfy the following conditions: 

1. To every pair of elements a and (5 in T, there correspond an element a + /3 in 
Tcalled the sum of a and fd, and an element a-(5 or a(5 in T, called the 
product of a and fd. 

2. Addition and multiplication are respectively commutative: For any a, (d in T, 

a + /d = /d + a a-(d = (d-a 

3. Addition and multiplication are respectively associative: For any a, (d, / in IF, 

{a + /d) + r = {a-/d)-r {a-/d)-r = a-{/d-r) 

4. Multiplication is distributive with respect to addition: For any a,P,y in T, 

a\ld^Y) = {a-ld)-v{a-Y) 

5. T contains an element, denoted by 0, and an element, denoted by 1, such that 
a + 0 = a, ha = a for every a in IF. 

6. To every a in IF, there is an element fd in T such that a + fd -Q. The element 
fd is called the additive inverse. 
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7. To every a in T which is not the element 0, there is an element y in T such 
that a-y = \. The element y is called the multiplicative inverse. 

Note 1: The symbol M and C is used to denote the field of real numbers and the field of 
complex numbers, respectively. 

Note 2: The set of integers and the set of polynomials do not form a field because they 
have no multiplicative inverse defined within the field. 


Definition D.2 (ring). More precisely: a commutative ring with multiplicative identity. 
A set with all properties of afield except for property 7 in Definition D.l. 


Note: The set of integers and the set of polynomials with real coefficients form a ring. 


Definition D.3 (linear space over a field). A linear space over afield T, denoted by 
{Ti,T), consists of a set, denoted by Ti, of elements called vectors, afield T, and two 

operations called vector addition and scalar multiplication. The two operations are 
defined over Ti and T such that they satisfy all the following conditions 

1. To every pair of vectors x, and x^ in Ti, there corresponds a vector x^ +X 2 in 
Ti, called the sum of x^ and X 2 . 

2. Addition is commutative: For any v,, Vj in Ti, x^ + X 2 = X 2 + x^. 

3. Addition is associative: For any v,, X 2 , x^ in Ti, (vj + ^ 2 ) +V 3 = Vj +(^2 + V 3 ) 

4. Ti contains a vector, denoted by 0, such that 0 + x = x for every x in Ti. 

The vector 0 is called the zero vector or the origin. 

5. To every x in Ti, there is a vector x in Ti, such that x + x = 0. 

6. To every a in T, and every x in Ti, there corresponds a vector ax in Ti 
called the scalar product of a and x. 
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7. Scalar multiplication is associative: For any a, (5 in T and any x in Ft, 
a{^/3x) = {^a/3)x. 

8. Scalar multiplication is distributive with respect to vector addition: For any a in 

T and any x^, in FC, a(^j+^ 2 ) = +a^ 2 - 

9. Scalar multiplication is distributive with respect to scalar addition: For any a, 
P in T and any x in FC, [a + P)x = ax + Px. 

10. For any x in Ft, \x = x, where 1 is the element I in T. 


Note 1: If J7'= R, then is called the n-dimensional real vector space. For 

brevity, the field of real numbers is implied and therefore is denoted as M". 

Note 2: [ p] is the set of all polynomials of degree less than n with real coefficients. 

Note 3: If a field is replaced by a ring M.., then (^FCMi) called a module over the ring. 
A nxl or Ixn polynomial vector can be considered as an element of the rational vector 
space , or an element of the module (^M"[p],R[p]^. A set of polynomial 

vectors is linearly independent over the field R(p) if and only if the set is linearly 
independent over the ring ]R[ p]. 


Definition D.4 (normal rank). The normal rank of a polynomial matrix 
A(p) e M""*" [p] is defines as e.g. 

1. The number 

rank A(p) = max rank A(p) 

2. The number of linearly independent columns (or rows) of A(^p) 

3. The number of invariant polynomials of the Smith form of A (^p^ 
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4. The dimension of the space spanned by the rows/columns 


Definition D.5 (rank of a rational matrix). The rank of a matrix M (^p^^ M”''" (ji?) is 
defined as e.g. 

1. The number of linearly independent rows/columns with rational coefficients 

2. The (normal) rank of N(^p) where N (^p^ is the numerator in a left or right MFD 
of M[p) 

3. The dimension of the space spanned by the rows/columns 

Definition D.6 (matrix pencil). A matrix pencil is a polynomial matrix M i^p^ of degree 
one. It is often written as 

^ {p) = pE + F 

where E and F are constant matrices and p is a differential/difference operator in the 
iMplace/Z domain in continuous/discrete time. 

Definition D.7 (regular matrix pencil). A matrix pencil M i^p') is said to be regular if 
M i^p^ is square and full (normal) rank. 


Definition D.8 (controllability index). The fact that the controllability matrix 

C = [5 ^5] 

has rank n means that there are n linearly independent columns in C . It could be that 
n such columns in the partial controllability matrix 

C^=[_B AB---A‘^-^b], l<q<n 
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can be found and that the smallest q is denoted as //. The controllability index ju of 
the pair is said to be the smallest q such that has rank n. 


Definition D.9 (unimodular matrix). A square polynomial matrix M(p) is called 
unimodular if its determinant is nonzero and independent of p. 


Definition D.IO (right/left divisor, right/left multiple). Consider the polynomial 
matrix equation 

A(p) = B(p)C(p) 

where A(p), B(p), and C(p) are polynomial matrices of appropriate orders. 

1. C(/)) a right divisor o/A(/)) 

2. B{p) is a left divisor of A(p) 

3. A(/)) a left multiple o/C(/)) 

4. A(p) is a right multiple of B(p) 


Definition D.ll (greatest common right divisor (gcrd)). A square polynomial matrix 
R(p) is a greatest common right divisor (gcrd) of N(p) and D{p) if R(p) is a common 
right divisor of N(p) and D(p) and is a left multiple of every common right divisor of 
N(p) and D(p). 


Definition D.12 (greatest common left divisor (geld)). A square polynomial matrix 
Qip) is called a greatest common left divisor (geld) of A(p) and B{p)if Q(p) is a 
common left divisor of A(p) and B(p) and Q(p) is a right multiple of every common 
left divisor Qfp) of A(p) and B(p). 
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Definition D.13 (right coprime). If a square polynomial matrix R(p) is a gcrd of 
N(p) and D(p) and R(p) is a unimodular matrix, then N(p) and D(p) are said to be 
right coprime. 


Definition D.14 (left coprime). If a square polynomial matrix Q(p) is a geld of A(p) 
and B{p)and Q(p) is a unimodular matrix, then A(p) and B(p) are said to be left 
coprime. 


Definition D.15 (Sylvester matrix). Let W(p) be a polynomial matrix of degree d such 
that 

W(p) = f^W,p‘ 

i=0 

with ^0 . Then the q*' order Sylvester matrix for W{p) is defined as 

ITo ITj ... W, 0 0 0 “ 

0 ITq 1Ti ••• IT^ 0 0 

0 0 0 
0 0 0 ITo IT; ••• 


S{W(p),q)^ 


Definition D.16 (characteristic polynomial). The characteristic polynomial of a proper 
rational matrix G(p) is defined as 

1. The least common denominator of all minors of G(p). The degree of the 
characteristic polynomial is defined as the McMillan or, simply, the degree of 
G(p) and is denoted by SG(p). 

2. Consider the proper rational matrix G(p) factored as 

G{p) = N{p)D-\p)^D-\p)N{p) 
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where N(p) and D(p) are right coprime, and N(p) and D(p) are left coprime. Then 
the characteristic polynomial of G{p) is defined as 

det D{p) or det D{p) 

and the degree of G{p) is defined as 

SG(p) = ^( det D{p)) = det D(p)^ 


Definition D. 17 (Right polynomial fraction). Every qxp proper rational matrix G(p) 
can be expressed as 

G(p) = N(p)D-fp) 

where N(p) and D(p) are qxp and pxp polynomial matrices, respectively. G(p) 
in this form is called a right polynomial fraction or, simply, right fraction. 


Definition D.18 (Left polynomial fraction). The dual expression of G(p) in Definition 
D.16 where the same proper rational matrix can be equally expressed as 

Gip) = D-\p)N{p) 

where D(p) and N(p) are qxqand qxp polynomial matrices, respectively. G(p) in 
this form is called a left polynomial fraction or, simply, left fraction. 


Definition D.19 (column/row reduced). A nonsingular polynomial matrix M(p) is 
column reduced if 

S (det M (p)) = sum of all column degrees 

It is row reduced if 

(5'(det M (p)) = sum of all row degrees 
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Definition D.20 (Bezout Equation). Consider two polynomial matrices A(p),B(p) and 
four polynomial matrices P(p),Q(p),R(p),S(p) such that 


Pip) 

Qip) 

Aip) 


~Gip) 

Rip) 

Sip)_ 

Bip)' 


0 


then the following can be said: 

1. A(p),B(p) are left coprime if and only if G(p) is unimodular 

2. R(p) and S(p) are right coprime 

3. B(p)Afp) = -S\p)R(p) 

This implies that if G(p) is unimodular, then right and left polynomial fractions of the 
same proper rational matrix are equivalent as 

H(p) = B(p)A\p) = A\p)B(p) 

where 

A(p) = -S(p) 

B(p) = R(p). 


Definition D.21 (Diophantine Equation). Given three polynomial matrices 
A(p),B(p),C(p) that are mutually right coprime and deriving equations in a Bezout 
Equation form, define the Diophantine Equation as 

C(p) = E(p)A(p) + E(p)B(p) 

from the top row of the Bezout equation and a left coprime representation such that 

R(p)A(p) + S(p)B(p) = 0 

from the bottom row of the Bezout equation. A solution Eq(p),Eq(p) exists such that all 
solutions (nonunique) can be expressed as 

E{p) = E,{p) + W{p)R{p) 

E(p) = E,(p) + W(p)S(p) 
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with W(p) an arbitrary matrix of appropriate dimensions. 


Definition D.22 (polynomial degree operator). Given a polynomial matrix M(p), 
define 

S^.M(p) = degree of eolumn of M(p) 

S^^M(p) = degree of row of M(p). 

Denote S^. the eolumn degree and S^. the row degree. 


Definition D.23 (right polynomial transfer matrix and control canonical form) 

Let a minimal realization of the system 


px(t) = Ox(0 + r u(t) (. 109) 

yit) = Cxit) + Duit) (-110) 

where x{t) e R"*, u{t) e R"", y{t) e R"^ and (0,r,C,D) are of appropriate dimensions. 
Beeause the system is a minimal realization, an equivalent representation of ( .109)- 
( .110) is the right polynomial transfer matrix form as G(p) = B(p)A (p). Let 
be the column degrees (i.e. the maximum degree of the polynomials in eaeh 
eolumn) of the matrix A(p) and p is the differential operator in eontinuous time or the 
time shift operator in diserete time. If z(t) is defined as the partial state [200] (or 
sometimes ealled the pseudo state [61]), then it is easy to see that it ean be related to the 
input u(t) and output y(t) as 


u(t) = A(p)z(t) 


(. 111 ) 


y(t) = B(p)z(t). 
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(. 112 ) 



The partial-state equations in ( .112) can be expressed in state space form by first 

defining the matrices 


i:(p) = diag\p^',p^\...,p gT?""""" 


(.113) 


L (p) = diag 


1 

P 


1 

P 

1^2- 


1 

P 


R 




(.114) 


and expressing the polynomial matrices of (. 111)-(. 112) as 

A(p) = 4S/p) + A,L^(p) (.115) 

B{p) = Bj:^{p) + B^L^{p) (.116) 


where are matrices of high degree coefficients whose dimensions are 

and n^.xn^, respectively, and are the matrices of all lower order coefficients 

with dimension x and n^xn^, respectively. The partial-state vector with dimension 
n, = p^+... + , can now be expressed as 


z(t)= (0,...(0 (.117) 


can now be expressed as 


z(t) = L^(p)z(t). 


(.118) 
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Substitute ( .115)-( .116) into ( .lll)-( .112) and assume is nonsingular to obtain the 
state space realization of the original system in (0.32) and (0.33) as 




A zit) + A^^uit) 


(.119) 


y(t) =(^i- ^hAh'A ) z{t) + B,Af^'u(t) 


(. 120 ) 


It can be shown [61] that the ensuing state space structure of (.119)-( .120) is in 
controllable canonical form for which the nonsingular matrix T can be found for MIMO 
systems using a method such as the one published in (Zak, 2003; [386]). 


Theorems 

Note: Theorems herein are stated without proof. Source information is provided. 


Theorem D.l ([61];Theorem G-6). A square polynomial matrix is unimodular if and 
only if its inverse is a polynomial matrix 


Theorem D.2 ([61]; Theorem 7.1). The controllable canonical form is observable if and 
only if D{p) and N(p) are coprime. 
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Theorem D.3 ([61]; Theorem 7.M2). A state equation (^A,B,C,D) is a minimal 
realization of a proper rational matrix G(p) if and only if the {A,B) is controllable and 
the {A,C^ is observable or if and only if 

dim(A) = (5'G(/)). 


Theorem D. 4 ([61]; Theorem 7.8). Let N(p) and D(p) be qxp and px p polynomial 
matrices, and let D(p) be column reduced. Then the rational matrix N(p)D \p) is 
proper (strictly proper) if and only if 

S^iN(p) < d^p(p) {S^.N(p) < d^p(p)) 

for i = 1 , 2 ,..., p . 


Theorem D.5 ([61]; Theorem 7.M3). All minimum realizations of G(p) are equivalent. 
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